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SYNOPSIS 


Sqn. Ldr. G.A.Y. Murthy, M.Tech., Aero. 
Indian Insl^i'tu'fca of Toch.no lo^y. Kanpun 
August, 1973 


PREDICTION OE BYNilMIC MTERAIi STiiBILITY 
OHilRilCTERISTICS 


Thv. equations of lateral-^Di recti onal motion of an 
airplane for the stick-fixed case are developed in terms 
of the mass, geometric and flight parameters. 


The coefficients in these equations are functions of 
of the lateral stability derivatives. Several methods of 
estimation of these derivatives and their values estimated 
by these methods are presented, J 


Computer programmes have been developed for estimating 
the time history of important airplane lateral motion 
parameters. 


j-i. sample case, based on the Cessna— 182 single— engined 
airplane, has been worked out in detail, covering the estima- 
tion of the stability derivatives and the prediction of the 
airplane motion in response to a disturbance, by solution 
of the lateral equations of motion. 



1. PTTROMQTIOH 


The calculation of the stick-fixed lateral dynamic 
characteristics of an airplane requires the solution of the 
lateral equations of motion. These are a set of three 
coupled, simultaneous, second order differential equations. 
The unknown parameters a,re, the angles of sideslip, roll 
and yaw and rates of roll and yaw, (Che coefficients in these 
equations contain stability derivatives and are functions 
of the airplane's mass, geometric and flight parameters. 

Two methods of solution of the equations of motion 
are considered. These are, the Laplace Transform method 
(Appendices A, B and C) and the fourth-order Runge-Kutta 
method (Appendices D and E), 

The Laplace Transform method has the following 
advant ages ; - 

(a) It provides the roots of the stability quartic. 

(b) The period and damping of various modes of 
lateral motion can be calculated. 

(c) It permits calculation of the free motion 
following any initial conditions, as well as, motion 
following application of external forces and moments, 

(d) It is very useful in further analysis such as 
frequency response studies by Bode plots etc. 



1. BTTROJUOIIOU 


The calculation of the stick-fixed lateral dynamic 
characteristics of an airplane requires the solution of the 
lateral equations of motion. These are a set of three 
coupled, simultaneous, second order differential equations. 
The unknown parameters are, the angles of sideslip, roll 
and yaw and rates of roll and yaw. The coefficients in these 
equations contain stability derivatives and are functions 
of the airplane's mass, geometric and flight parameters, 

!Ewo methods of solution of the equations of motion 
are considered. These are, the Laplace Transform method 
(Appendices A, B and G) and the fourth-order Runge-Kutta 
method (Appendices D and E), 

The Laplace Transform method has the following 
advantages:- 

(a) It provides the roots of the stability quartic. 

(b) The period and damping of various modes of 
lateral motion can be calculated, 

(c) It permits calculation of the free motion 
following any initial conditions, as well as, motion 
following application of external forces and moments, 

(d) It is ve 2 ?y useful in further analysis such as 
frequency response studies by Bode plots etc. 
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Solution obtained by Runge-Kutta method provides 
only the time history of the unknown parameters. This method 
Was considered for its simplicity and for comparison v/ith 
the results of the Laplace Transform method. 

Several methods of estimation of lateral stability 
derivatives are available. These are presented and discussed. 

It Was considered very important to take up a sample 
case and make a complete analysis, Oessna-182 (Pig, 1 (a) , (b) 
and (c)) Was chosen for this purpose because of the availa- 
bility of this model in the Depart;mont and the possibility 
of someone else pursuing this programme by obtaining flight 
data for correlation with the data estimated in the report. 

The motion of the airplane as calculated here is in 
the form of time history of angles of sideslip, roll and 
yaw and rates of roll and yaw following any initial conditions. 

The results CEin, by proper modification of the pro- 
gramme, be extended to calculate motion resulting from impr- 
essed forces and moments which are functions of time, 

/ill the forces and moments in the text are referred 
to stability system of axes which are defined in Pigg. 2 and 3. 

Ts^ical input /output for the solution of lateral 
equations of motion and prediction of time history of angles 
of sideslip, roll and yaw and rates of roll and yaw are 



2» _5QU'ATI0]^S OP L /i(PTi}R/.T, 

IpTIOHS Qg AH AlRPT./iTTO 


2.1 The dimensional form of the equations of lateral 
motion (stick-fixed) arej 

( 2 . 1 ) 
( 2 . 2 ) 

’ 8? § - if - H § - &ift)(ten',‘ )| 


m 


2. 

dt 


^2 at 2 
2* 


3I 

9r 


dt 5T ^ 


0 




^ ^ aiT 

3r df ~ 


0 


+ 




Y 


c 


0 


(2.3) 


Dividing equations (2.1) and (2.2) Dy ^ P y^s^b and 

(2.3) by ^PY S, the non dimensional form of equations are 
obtained. 





at + 

ds 



1 






4 


irt 

where s = non-dimensional time 

Equations (2#4) are multiplied by (iV/'Sb), the rej.a- 

tive density factor ji and are written in the following form 

(D^-l^D) ({) + (K^D^-lj,])) f - lpP-1^ = 0 

(KgD^-n^DT (|) + f - npP-n^ = 0 / (2.5) 

(-ypD-Cj^/2) + (D-y^D~aj./2 (tan^)) if + (l>-yp) P~y^ = 0 

where D = d/d<r' and the other coefficients are explained 
in the list of symbols, 

2,2 Solution of Equations of lateral Motion by Laplace 
Transform MethodJ- 

TaEing Laplace transform of equations (2,5)j we have, 
with the given initial conditions, 

- ipA^ ipAp^ = 

If^ - npAP;s ^ ^2 

' 

^ _ r 

(-y^A^ - ^X) Ij)^ + |a^ - y^X^ - Cj./2Ct£tn’<)A| 

% 

+ (A^ - 7pA ) = r. 
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Where r, = (A ^ ^ ) lf„ 

+ A(Di|))^ + K^A Cuf)^ + 1^ 

^ 2 ~ (KgX^ - n^) (|)jj + ( - rij, A ) 

+ K2A(Df))^+A(]l|r) +n 

W O 

Tj = - y-pA + (A - y^A ) + A Po - 

and L ((j)) i etc. 


Solving equations (2.6) {ip^ , and P;, can be 
expressed as 


A _ ^qA + a.|A -f a2 A -f + ^A + a 

~ 3 A* -v 4 Ti 2 — ' 

A (AA^ + By + 0 A + dA+ e) 


= ^0^ + b.|A + bgA ^ + b^A ^ + b^ A+ bp 


A (aA + bA^ + gA^ + dA+ e) 


'A 


+ C..A'^+ C.,A‘^ + C' 5 :A+ q. 


o + 


^ 2 
+ Cg 


A(aA^ + bA^ + gA^ + 1)^4- 


E) 


( 2 , 8 ) 


(2.9) 


The coefficients a^, a^, c^,..o, and 


5' 

A, B, G, D and E which are the coefficients of the 
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characteristic equation, 

A'^Vt + bX^ + gA^ + -D/\4- E 


0 


(2,10) 

\ 1; ’XA 'I*'"'- ■' ''' 


are functions of airplane mass geometric and aerodsmamic 
stability derivatives, and the expressions for these coeffi- 
cients appear on pp,^r^3( Appendix A), 


!Ehe expressions for l|r and p can now be obtained 
by taking the inverse Laplace transform of equations (2,7) 
to (2,9). 


Assuming that all roots 7\j_ of fbe characteristic 


equation are distinct, the inverse laplace transform gives 

A -i 

{|) = A^e 1 + A 


2 ^ A 


+ A^e 


4- ^ 4. 4 “ Ag 


( 2 . 11 ) 


o-'A. o-^/i o-A. 

f = B^e ' + B^e ^ + B^e ^ + B^e t ^ Bg(2.l2) 

cr^A. 0-"Xp cr^A o'^A. 

p = C^e + G^e ^ -f G^e ^ + G^e ^ + G^ (2.13) 


and the expressions for rate of roll, p and rate of yaw, 
r, are obtained by differentiating equations (2,11) and 
( 2 , 12 ) respectively, 

i bX . cr>. 

P = (A^A^e + ^ 2 ^ 2 ^ ^ ^ ^^ 5 ) 


( 2 * 14 ) 
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^ + B^h^e ^ + b^>\^ e 


^ + B^A ^ e° \ b ^) 


( 2 , 15 ) 


03ie expressions for the coefficients A-j,.,Ag, B^,.,,Bg, 
appear on pp. ^4’(Appendix A). 



3. ESTIMATION OP TIME HISTORY OF LATERikL MOTION 


3.1 The process of obtaining the lateral motions of an 
airplane can be divided into four stages. 


Stage l i The values of the following parameters are to be 
det ermined/estimat ed j 


(a) 

(b) 

(c) 


Mass characteristics; m, k , k , n and p 

\ i ■ 

G-eometric characteristics; Area S, span b etc. 
Fli^t conditions; velocity ?, lift coefficient 
and fli^t path angle . 


(d) Aerodynamic Stability Derivatives. 

i) Side-slip derivatives On > 0 , Oy 

■^p p ^3 

ii) Roll-rate derivatives C-, j 0 , Gy 

P P 

iii) Yaw-rate derivatives C-, , 0 , Cy 

*^ 1 ' 

The methods by v;hich these derivatives are estimated, 
are elaborated in the next section, 

(e) In cases where impressed forces are used as dis' 

turbances, the values of factors On , C , Gy appropriate 

c ^c c 

to the particular problem are to be determined. 


Stage II ; Prom the above known quantities, the non-dimensional 
parameters as used in the calculation of motion are evaluated. 
Prom these parameters, the coefficients A, B, 0, D and E 


^ag , e III ; The roots of the above characteristic equation 
are then obtained. In the programme used, the method of 
reduction to cubic has been used, 

§:faage IJi The appropriate factors 0^,,.. 

are evaluated. Some of these factors may be complex conjugat 
pairs depending on the nature of the roots /A. of the charac— 

JU 

teristic equation, Accordingly the equations of motion are 
written as follows; 


(a) All Roots Real 


E A^e + Ar-cr' + a 


1=1 


t = S B . 

• -L. 


i=1 


cAA. 

0 ^ + B^o--' + Bg 




4 cr'A, 
1 


^ S 0. e + G 


i=1 


P = 


= X 


5 

S A.X . e ^ + Ar 


i=1 


1 1 




4 . 

V 1® + 


(3.1) 


(b) Ono Gomplox Pair of Roots and Two Real Roots 

At, = R + II j 

f 

A2 ^ R - il j 


complex conjugate pair 
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K_^e 0os(crl + w^) + 5 + A^e + A^cr + A^ 

a->_ cfA 

V 008(0-1 + iVj) + BjO ^ * B e + B O- + B 


= K e Oos{cKl + w_) + 0. 


vua^,o X + w^; + o^e + C^e ^ 


(3.2) 


2 


~ K /R +I e Oos(cr'l + w + tan z) 
^ A R 


o-'A o-^ 

+ A^A ^ + A^A^e 


1 I / 2 2 -i-l t 

r ’*’^ ® OosCcf'I 4- w^ + tan i) 


o'X o-'X -| 

■"Vs® -"^'4® ■" ^5 


Vi? here 




2x/R‘r + I 


-1 J- 

tan 

""a 




(3.3) 


Kq = 2 Aq + 


ri 52 


Ai = + I^i 


i A, = R, - I.l 


■A "A" J 


The steps to determine R^ and have been shown in 
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(c) Two Pairs of Oomplex Itoots; 

R + il 

Aj’ ^4 " = t il' 


cr'Ki 


V Cos (tr'I + wp + K^e Cos + *•) + 


A' 


o^llt 


Kg Oos (cr^I + Wg) + K^e Oos (cr'I* + w») + Becr^+ 


Kq© Gos (o'' I + Wq) + K^e 


- o-'iii 


Oos (cr'I* + w«) + Cj 


1 / 2 2 ^ -1 
^ r ® ^os (o-l + w + tan |) 


g; + 


+ IC^/Rt +It^ e 


e Cos {cRI» + w* + tan ^.) f 

A gt/ 


.2 2 


^ V K +I e Gos (cr" I + + tan |) + 


/-2 <^R* - -•) - 

+ K^vAi* +lt e Oos (cKl* H- w| + tan £1) | 


' ' ' (5.4) 

Where quantities with primes are defined in exactly the 
same way as in (3.3). 
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The appropriate equations of motion (3.1 ), (3.2) or 
(3.4) can now bo solved by substituting the values of non- 
dimensional time factor cr^ to obtain the time history 
of lateral motion parameters, p, p or r, 

3 . 2 Motions Resulting from Arbitrary Disturbances t 

The motions resulting from arbitrary forcing functions 
can be obtained from the motions resulting from constant 
impressed forces and moments by methods explained in Ref .4. 

In Ref. 4 it is shown that although component notions 
of an airplane must be calculated simultaneously, the effects 
of component disturbances may, by the principle of super- 
position, be calculated separately and later added in desired 
proportions. This principle of supeiposition is elaborated 
with examples in the same reference. 



4> ESTIMiglOM OE LATERAL STABILITY DBRIYATIVES 


G-EMERAL ; Methods of estimation of contribution of all 
the main components of an airplane to all the ' derivatives, 
for aiii)lanes having any sweep angle or aspect ratio, have 
been presented in various publications. 

As long as these methods are based on available flight 
test or wind tunnel test data on similar designs or theo- 
retical or empirical formulae, the values obtained are suit- 
able only for first approximations of dynamic stability. 
Extensive wind tunnel test data on scaled or full scale 
model of the particular configuration will be essential for 
accurate estimations. 

The estima.tions made here depend on theoretical treat- 
ments, based on assuinptions of steady roll and yaw or experi- 
mental data obtained principally from tests made under 
conditions of steady roll and yaw* 

In the present study, the effect of Mach number on 
the derivatives is indicated for the subsonic case. But 
for the sample case of estimation of stability derivatives 
for the Gessna-182 airplane, the effect of Mach number is . 
neglected* The effects of power are briefly treated in Ref, 4, 
but these effects are also neglected in the present study. 
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The estimation of stability derivatives has been based 
on the available formulae and methods of Refs. 1 to 5* 

The first three are text books, which have drawn material 
from published and unpublished research papers. The fourth 
is a summary of the available methods of estimation. The 
fifth is a veiy recent publication on the subject and has 
drawn material generously from ’USAR Stability and Control 
Datcom* v/hich is one of the most reliable sources in use 
in IJ.S.A, for estimation of stability derivatives. 

The formulae and methods contained in the above five 
references have much in common, but in Ref. 5 additional 
(empirical) correction factors are provided and also geome- 
tric twist is accounted..! or. Wherever there are wide varia- 
tions among the formulae, it has been so pointed out and 
the results obtained by these various methods are tabulated 
in Table 6 (a) and (b) for the two fli^t conditions consi- 
dered, The values of derivatives which are provided by the 
manufactures of Cessna airplane are also given in the same 
table for comparison. A discussion of these results follows 
in section 5. 

4 . 2 The Side-slip derivatives Q y > 0^^ , s - 

p p p 

These derivatives may vary widely and often in an 
unpredictable manner with angles of attack and hence no 
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essential to base these estimates on wind tunnel test data. 

In the absence of such data, reasonably accurate estimates 
Can be made by applying suitable corrections, to available 
data for similar designs. Certain theoretical and empirical 
formulae have been evolved for estimating the contribution 
of the major airplane components. 

The side force derivatives are affected by such compo- 
nents as lea.ding edge hi^-lift devices, stall control devices, 
trailing edge flaps, nacelles, external stores like fuel 
tanks etc., canopies and dorsal and ventral fins. Generally, 
v/ind tunnel test data is essential for estimating all these 
effects. Several research papers dealing with the effects 
of these individual components are available for example 
in Ref ,4. 


4.2.1 . Side gorce herivative Due to Side-slip , Cy 

P 

4. 2. 1.1 This derivative gives the variation of side force 
coefficient with a positive side-slip angle p. Cy is usually 

P 

negative. This derivative contributes to the damping of the 
Dutch Roll mode; thus, large negative values of Gy might 

P 

seen desirable. But large negative values of Gy , however, 

V P 

may create a large time lag in the airplane’s response and 
cause it to react sluggishly to the pilot’s commands. The 
main contributions are from the body and the vertical tail. 
For v/ings with dihedral, the contribution of the wing may 
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compared to that of the vertical, tail* The body contribution. 

arises mainly because of separation of flow over the fuselage. 

The main contribution to 0^ is from the vertical tail. As 

P 

far as possible an accurate estimate of this contribution 
is essential, since this factor is used in estimation of 
-vertical tail contribution to several other derivatives. It 
is specially significant at low angles of attack. It is 
highly desirable to have tail-off and tail-on wind-tunnel 
test data for an accurate estimation* 

. Wing plan form and wing location can cause large diff- 
erences in side-vi/ash and dynamic pressure at the tail. It is 
difficult to predict the side-wash effects accurately and 
hence wind-tunnel test data is very helpful, 

4. 2.1.2* The contribution of various components are estimated 
as follows: 


t Cly (4.1) 

P Pw Pb h 

(a) Wing: Generally thi a. contribution is small. For wings 
with large dihedral, this contribution is estimated from 
the empirical formula given in Ref. 5. 

0 =- 0.0001 |f |.57.5 Crad“'') (4.2) 

Pw 
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wiiere 

I != Geometric dihedral, deg, 

(b) Body t- A rough estimation of this contribution as 
given in Ref. 2 is 

(V = - 2 

Pb 

where S-j^ is an equivalent base area and S is the wing area. 

lies betvi?een zero and (%/4-)iL where h is the maximum 
fuselage depth. But it is very difficult to obtain an exact 
, But since 0 is small compared to vertical 

Ip 

2B 

tail contribution, a high degree of accuracy is not needed* 
'is independent of Mach number as per slender body 
the^y, 

A better approximation for this contribution is given 
in Ref, 5 as follows 


value of 




Gy = - 2K^ (Sq/S) (rad"^) (4.4) 

Pb 


where 

= a wing-body interference factor obtained as a 
function of the wing location with respect to 
body i»e., /(d/2), from Pig, 4. 

where . 

Co 

s= the distance from body-centre -line ..quarter chn-ra 
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the quarter chord point helow the body centre 
line* 

d = naximum body height at wing-body intersection, 

= cross-sectional area of the fuselage a,t the 
point along the body where the flow ceases 

to be potential, 2he distance X is a function 

o 

of X^ , the body station where dS^/dX first reaches 

its maximum negative value. She distances X and 

o 

X^ are related as in l^ig, 6, To find and X^, 
variation of fuselage cross-sectional area along 
the body centre line has been tabulated in Table 2, 
and plotted as a function of X.. (Fig,5), X^ is 
determined from this figure and corresponding value 
of from Fig, 6. The value of is then read 
from Fig, 5* 


(c) Vertical Tail i A theoretical estimate of this contri- 
bution is made from the formula given in Ref ,5* 



kO^ (1 + d(rydp)fj^ S^/S (rad"^) 

tty 


( 4 . 5 ) 


where: 

k = an empirical factor depending on the constant (by/2r^) 
(defined in Fig, 7) andTly and is obtained from Fig, 7, 

0^ = tail lift-curve-slope based on the effective aspect 

ay 

ratio which is approximately 1,55 times the geometric 
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aspect ratio (b^/s^) (Ref.3), and is calculated as 

follows 


27i;/i 


"7 


eff 


a 


7 


2 +'/ 


/ 4 P' 


(4.6) 


eff 


Jl 2 


-0 . v^- 


P' 


^) + 4 


where! 

P = /1 - 

= ratio of actual average tail section lift-curve-slope 
to 2% 



K 



(4.7) 


^atio of the aspect ratio of the vertical panel 
in the presence of the body to that of the isolated 
panel and is obtained from Pig, 8 as a function 
of "I ^3i(i /\ y • 


“^(B) 


ratio of the aspect ratio of the vertical panel 

in the presence of the horizontal tail and body 
to that of the panel in the presence of the body 
alone, This ratio is obtained from Pig, 9 as a 
function of Z/d, and rt 
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defined in the sane figizre) 

Kg = a factor accounting for the relative size of the 

horizontal and vertical tails and is obtained from 
Mg* 10 as a function of Sg/S^. 


The effect of sidewash and variation of dynamic pressure 
at the vertical tail is calculated from the empirical 
formula 


(1 + doVdp)’';^ = 0.724 + 3.06 


1+OosAq^^ 


0.4 



.009 A 


(4.8) 


4.2, 1,3 References 1 to 3 suggest utilisation of wind 
tunnel test data for estimation of C„ • According to Ref. 4 

P 

interference effects are so large that a generalized formula 
would not be completely satisfactory. A method of applying 
correction to the data on a similar design ha.s, instead, 
been recommended, Glie method given above from Ref, 5 probably 
gives the most accurate values since interference effects 
based on experimental results are included. 


4,2.2 Variation of Riolling Moment Coefficient y/lth 
Sideslip Angle , 0, 

l.i - 

4 .2, 2.1; This derivative is also called the * dihedral 

effect*. C-, may be positive or negative and for some 
in 
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configurations it has large variations. It is a very 
important design parameter. It aids in damping both the 
Dutch Roll mode and the spiral mode. Small negative values 
of 0-| give favourable Dutch Roll damping characteristics, 

P 

but for spiral stability large nega,tive values are necessary. 

The main' contributions to 0^ are from the wing-body, 
horizontal tail and the vertical ^ail. 

(a) Wing- Body Contribution ? This contribution is again due 
to three aerodynamic effects, 

i) Wing geometric dihedral 

ii) Yfing sweep angle 

iii) Wing position on the fuselage. 

(i) When an airplane with positive dihedral side-slips, 
the right wing will be at a liigher angle of attack 

(/la s= p/”") than the left wing. The increased lift on the 
ri^t wing and the corresponding decreased lift on the 
left wing, cause a negative rolling moment. This rolling 
moment is proportional to the geometric dihedral. 

(ii) Y/ing sweep angle has the same effect a,s geometric 
dihedral. This effect is proportional to the lift coefficient. 
For large i.e., at low speeds, 0^ duo to sweep can be 

P 

very large negatively and becomes small negatively at high 
speeds. This poses a serious design problem for hi^ speed 
airplanes with sv^ept back wings. 
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iii) The fuselage contribution to 0-, is also very signi- 
f leant. The flow field is considerably modified by the 
shape of the body. The cross flow component of the stream 
induces vertical velocities which , combined with the main 
stream velocity alters the local angle of attack of tho 
wing. Consequently, a hi^ wing configuration produces a 
negative rolling moment and vice versa. The effect also 
depends on the fuselage length ahead of tho wing, its cross- 
sectional shape and planform and location of the wing, 

(^) Horizontal Tall Contri bution, The horizontal tail 
contribution can be explained in exactly the sam e manner 
as for the wing. Except when the horizontal tail size is 
considerable, this contribution is generally negligible, 

(c ) V ertical Tail Co ntribution, The sideslipping airplane 
produces side forces which give rise to significant rolling 
moment if the vertical tail aerodynamic centre is offset, 
considerably from the rolling axis (i.e,, if z is large). 

4*2'»2*2 i The estimation of the contribution to G-, of 

p ' 

various airplane components is carried out a.s follows; 

C-i = G*i t 0- ( 4 . 9 } 

P - Pl,^p pjj Py 

(a) Wing-Bo dy; The wing-body contribution is calculated 
by the following formula (Ref ,5). 
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Cl 57.3 




“if ^ 

X A I 


J ^lo ^ ^1 'I 

+ I i(y=^) % + -jJ-l + (A% ) 

■ r ' J a 


w 


A Cl 

+ 0 tan A — P 

0/4 0 tanAo/4 


rad 


-1 


(4.10) 


where: 


\ b '^ ^X airplane steady state lift coefficient 

effect of wing sweep contribution and is 
obtained from Hg, 11, 



- compressibility correction to sweep and is 
obtained from I'ig, 12, 

= a fuselage correction factor obtained from 
Pig. 13. 

“ contribution due to aspect ratio and is 
obtained from. Pig. 14 , 

= wing dihedral effect obtained from Big, 15, 

compressibility correction factor to dihedral 
and is obtained from Pig, 16 , 
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body induced effect on the wing height as a 
function of the cross-section of the body 
and is given by the empirical formula. 






r 

where 

d 

av 


-0.0005 A (-f^) deg. 


(4.11) 


/ 4(average cross-sectional area) 



(4.12) 


The average cross-sectional area is calculated from Table 2 
as a weighted average. 


another body induced effect on the wing hei^t 
as a function of a and d _ i.e., the location 
of the wing w.r.t. fuselage and is 'given by 

(ACi^) 2 ^ = (deg"*'’) (4.13) 

AO 

■^G • ' . 

( = a wing twist correction factor obtained 

tanA^/^' 

from Fig. 17 . 

Refs. 1, 2 and 4 give a formula for the wing-body 
contribution which is almost similar to the above. 


(‘AOi ) 

%i 
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0. 


(h. + W ) 

w ^ a V av 


\ = “1+ ° i - r * 1.2 ir^ 

P\ 7 B I* P 


/ 

0. 


(4.14) 

where the values of (c^)> Oj are given in graphs and 

L pp 


^av fuselage height at wing root 


w 


av 


average fuselage width at wing root. 


w 


The formula of Ref. 5 takes into account the effect of -sweep 
and geometric twist* 

Ref • 3 gives the following formula for the wing-body 
contribution 




0 


(4,15) 


The wing dihedral contribution (O^ is given by 

P 

^'ap5w= (4.16) 


(A.0, ) . is the contribution to 0-, due to wing-fuselage 

P ^ -^P 

interference and is given below, 

(<^(^) 

Hlgb wing -.0006 

Mid wing 0 

Low wing + . 0008 
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Is the contribution to 0;^^ due to wing-vertical 
tail interference and is given below 




High Wing 

.00016 

Mid wing 

0 

Low w ing 

-.00016 

4 '*■,)= 0 1 _ a;, 

+ .0^1 - (rad”'’) 


“ (4.17) 

where 



k - 1,0 for straight wing tips 

k = 1.5 for roimd wing tips. 

Kie above formiila is a revised one for unswept, 
elliptical Wings with zero dihedral taking into consideration 
changes in taper ratio, for the wing contribution to 0 . 

h 

(b) Horlsonta_Iail: The contribution of the horizontal 

tail is Calculated from 


h I 


q >) 

rn H 

' ST“ 


') rad 


(4.18) 


WhPrP O n n 

Ip ® calculated by treating the fuselage- 

horizontal tail in the same manner as fuselage-wing has 

been treated. Significance nf r* j n 

. ignriicance oi C depends on the magnitude 

of (S„ VSb). 
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(°) Eerticai^: Ihe vsrtioal tall contribution io 
estimated from the following foimula (Ref, 5 ) 


% 

Py 


55 

Mm 

b 


(z Cos 


a 


= Gv • - ~ ^ 

Pv 


(4.19) 


where and are as defined in Pig, le. 


Refa. 1 and 5 give the following sli^tly different 
forimila for G. 


'1 


Pi 


2 

~ •'Iv 


(4.20) 


Of Ta wing Moment Ooeffloient w,-en 
Sideslip Angle , c . 

np 

4.2,3. 1 This Is also oallod the 'statio directional- or 
eather cook' derivative. It is a measure of the tendency 
of the airplane to align itself in sideslip, like a weather 
cock, with the relative wind. It Is one of the most impor- 
tant parsneters in design. The value of 0 ^ determines 
primarily the I^toh Roll natural frequenoyPand affects the 
spiral stability of the airplane. It varies considerably 
with tall size. For stability and good flying qualities 



28 


0 should he as high as practically possible. Ihe vertical 
(3 

tail, fusela,ge and wing contribute to G , but the tail 

Up 

contribution dominates. Por positive sideslip, the vertical 
tail causes a positive yawing moment and the fuselage 
contribution is generally a negative yawing moment. The 
wing contribution is usually positive but generally .neglec- 
ted, The value of 0 must be positive for stability, 

4. 2,3, 2 The contribution of the various components to 
G is estimated as follows 





( 4 , 21 ) 


(a) Wing ; Wing contribution is generally very small. and 
neglected. However, Refs, 1 and 4 give the following 
foimula for 0^ 

h 

'w ^ 

- ^n- 

whereas. Ref, 2 gives an additional term 

0 . . 0 .„ 

where, (Q-p )” 0,075 and ("-^) can be obtained from 


Ref. 4 



3. gives the following formula 


0 . 


n, 


'W 


A 1/2 

= .00006 (A) deg"'' 


( A in (A eg j 
Which for unswept wings is zero. 


(4.23) 


(b) Bo^i Ihe fuselage contribution based on slender body 
theory gives (Ref, 2) 


G 

Up 

^•’Pus 



(4.24) 


where ?ol^ la the volume of an equivalent circular body 
having the same side-view as the real fuselage, i.e., 


Volf = ^ / d^ dx X I 3 (4.25) 

But the slender body theory tends to oyer -estimate 
the fuselage instability, since it neglects the eeparation 
at the rear of the body. The body contribution including 
the interference effect of the wing on the body is ostima- 
ted from the formula (Ref. 5) 


Onp^= - 57.3 i (rad-l) 

(4.26) 

where 

Kj[ = an empirical factor for body and body and 

Wing interference effects and is obtained from 



Fig, 19, The factor depends on the ahape 
of the body, position of the 0,G, , the cross- 
section of the body etc. All these factors 
for Cessna are shown in Table 3, 

= a Reynold’s number factor for the fuselage 
and is obtained from Fig, 20. 

- body side area 
S 

Ig = length of the body 

The body contribution given above is valid for all subsonic 
Mach numbers. 

Refs. 1, 2, and 4 give the following formula for 
the body contribution 





where 


3t 



an empirical constant obtained from fig* 21 
as a function of and 


^CG 


= the distance from nose to C.G, 
lig# are defined in fig. 21. 


(c) Vertical fail : The contribution of vertical tail to 
0 is estimated as follov?s {Ref. 5) 


c 


n. 




(l^ Oos a + Sin a) 
b 


(4.29) 


4»3 Roll Rate Derivatives : 

When an airplane rolls with angular velocity p, 

about its x-axis, the motion of its wing tip is helical 

pb 

?/ith a non-dimensional helix angle equal to ~ . This 
motion affects the airflow along the wing and tail surfaces 
as illustrated in fig*22, for two points, viz., wing tip 
and fin tip. These angle-of-attack changes introduce 
perturbation in the a,erodynamic forces and mom.Gnts, 
further, the changes in the wing load distribution causes 
a modification of trailing vortex sheet. The vorticity 
distribution is no longer symmetrical and a sidewash is 



indue Gd on the vertical tail. Further, the trailing vortex 
sheet is no longer flat but helical. These sidewash effects 
ha,ve been studied theoretically and experimentally. 


Fa-ri n_ of Side Force Coeff icient with 

Roll Rate, Cy ; 

P 

4. 3. 1.1 This is not an important derivative in estimating 

the dynamic stability cha,r act eristics of an airplane. The 

main contribution is from the vertical tail, even thoufii for 

some configurations the wing contribution may be significant. 

The si.G:n of C may bo either positive or negative. 

P 

4. 3.1. 2 Tlie main contribution to is from the vertical 

P 

tail and is estimated from the rel,ation 


Cy O-G 

■ P Py 


? “ G 




(z Gos a - 1 Sin a) 

= 2 Gy^ (4.30) 

Py 


In the above relation, the effect of side?/ash is 
neglected. Ref. 4 gives the following modified relation 
which approximately accounts for the effects of v/ing side- wash 



2 
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Cos a - 1^ Sin a - 
_ _ 


■'V 


c 


Y, 


Py 


(4.31) 


If the effects of sidewash can be estimated the following 
relation from Ref. 1 can be used. 


^Y 


0, 




Py 


"a. 


Y 


L- 


I - Cos a - 1^^ Sin a) + 


9 c ■' 


V 


V 




(4.32) 


All estimation of the effect of sidewash can be 
made from the graph in Ref. 6. 

Refs. 2 and 3 neglect the value of G„ » 

Xp 

4,3.2 Ya-riat io n of Rolling Moment Coefficient with Roll 
Rat e j Ct 

■^p 

4, 3. 2.1 This derivative is also known as the ‘ dojnping“in~ 
roll’ derivative, since, it represents the resistance of 
the airplane tc rolling. When an airplane rolls, the angle 
of attack on the down-going wing is increased and that on the 
up-going wing is decreased. This changes the lift distribution 
on the two wings and a rolling moment, which opposes the 
initial disturbance, is produced. 

Cn is a negative constant so long as the local 
■^p 

angle of attccck remains below the stalling angle. If it 

exceeds, C-. may become zero or even positive and the wing 
P 

may 'auto rotate’. 
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The main contribution to C-, is from the wings. 

P 


In case the horizontal tail is very la,rge, its contribution 
may not be negligible. The contribution of the vertical tail 
is usually sma,ll. 0 is the principal determinant of the 


damping-in-roll characteristics of the aircraft* 


4.3. 2. 2 The contribution to 0^ of the various components 

P 

is estimated as follows 


= Cl + Gi + Gi (4.33) 

P PvifB Ph % 


(a) Wing-Body ; The wing -body contribution is estimated from 
the following relation (Ref, 5) 


C 


1 




(4.34) 


where 


P Ot 
(— ^) 


roll damping parameter obtained from Pig, 23 

-1 

as a fimction of A = tan (tan /\^^/p) (deg. ) 



ra,tio of average wing section lift curve slope, 
C-j , to 2tt; . 



Based on strip-theory t echni qu Ref , 2 gives the 
following relation for the wing contribution 


0-. 




■R 


12 


Yl 


1+ /\ 


(4.35) 


Shis generally provides a higher negcative value for 


■p 


G 


1 


= (^1 ) 

Gj^-0 (Gt )q Q 


Another relation given in Ref. 4 is as follows 
'a)G_ 


(%), 




^ o2 i 1.2 Sin^A I 

^ ^ I (A+4 CosA)( 


1. C G — -J* ) 

8 7CA ^ 


(4.36) 


where 


(R 4 =0 

p V 


<%4^=0;a=2* 


'"‘_A+4 Cos /\ 
A 

L. X 


■j + 4 Co s / \ 


(4.37) 


and (G^ )q _q. g_2^ is obtained from charts of Ref .4 as a 
function of aspect ratio, sweep angle and tamper ra,tio. 

Ref ,3 reproduces a chart from NAGA-TR 635, which 

gives the values of On as a function of aspect ratio and 

% 

taper ratio, 

(^) Horizontal Tail : She contribution to On of an extremely 

P 

large horizontal tail can be significant. It is calculated 
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from the relation (R’ef ,5) 


Stt 2 

“l = 0-5 ChjH <r) 


Ph 


P 


(4.38) 


where (0^^ )pj is the valLie of G based, on the horizontal tail 
P -^p 

geometry and is found from equation 4.34. The factor 0.5 is 
introduced to account for the rotation of flow produced by 
the wing (Ref, 4). 

(^) yertical Tail ; The contribution of the vertical tail is 
estimated from the relation (Ref. 5) 


0 


1 


P' 


V 


z-n- 2 

2 (^) Qj 


V 


(4.39) 


The above foimula is modified (Ref .4) to allow for approximate 
correction for the effect of the wing on the dajnping~in~roll 
of convent ionali/lo cat ed vertical tail surface 



2(z^ Cos a - 1-^ Sin a) 

- — — — — — (zy Cbsa~l,yffl.na~z,^) 


By 


(4.40) 
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4.3.3 Variation of Yawing Moment Goefficient with 

Rol l Rate , G i 
P 

4»3*3k1 This is one of the cross derivatives that cause 
coupling between the rolling and yawing motions of an airplane. 

The main contributions are from wing and vertical tail. The 
wing contribution is a negative yawing moment. Tie vertical 
tail contribution may be positive or negative depending on 
the tail geometry, the angle of attack and sidewash from the 
wing, C affects the Dutch Roll damping in that the larger 

Up 

its negative value, the less Dutch Roll damping. Thus, a positive 
C is desirable. 

The wing contribution consists of two parts. The first 
part is due to change in profile drag due. to change in wing angle of 
attack^ Dor a positive roll, the angle of attack increases on 
the ri^t wing and decreases on the left wing, which normally 
causes an increase in the profile drag on the former and a 
decrease on the latter, thus producing a yawing moment. 

The second part is due to the fore and aft inclination 
of the lift vector caused by rolling. The lift vector on the 
rl gbt wing is inclined forv^/ard and that on the other, backward. 
This causes a negative yawing moment proportional to 0^^* 

4, 3. 3,2. The contribution of wing and vertical tail to G 

P 


is estimated as follows. 
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'11 = + C 




(4.41) 


(s-) The wing contribut 


ion is ii:t iwo pctnij 


(AC^)^ 




(4.42) 


= (C ) + (G ) 

^ 1 ^ n ^2 


inip 'theony using tv\'o diinGnsiortpii -fi -nv 4 . • 

-A^uensional flow assumptions lead to 

the formula (Ref, 2) 


(/_»^)i 


( Ot, ) 


-L- (ItLd.) 

12 S+ A ^ 


(4.43) 


( Ao ^)2 


(\)a 


L ) 

12 (1+A ) 


(4.44) 


These formulas are good only for a first approximation 


ana give an over estimation of (o„ even for unwept wings, 
low aspect ratio and high sweep angles, the effect is 


more pronounced. The estimate for (0 >2 is satisfactory 


except for low aspect ratio wings. 


Ref. 4 provides emplrioal charts which give a more 

ac curat e ent i mpit. i nn )l (^C„ ). 
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from Pig. 24. 


Ref. 5 provides the following relation for the wing 
contribution. This formula takes into account the wing sweep, 
and geometric twist. 



G-j tan a 
% 



C 

n 

tan a - ( 77 -^) 


0 t -=0 * 

if 



1 




n 

a 

o „ f 


f 


6 

f 


(4.45) 


where J 


C 


1 


% 


wing contribution to 
(4.34) 


C-j calculated from 
P 


0 

M 


wing angle of attack 
wing lift coefficient 

total airplane lift coefficient. 

variation of the yawing moment coefficient 
due to roll rate at zero lift, given by 


«n 


Cos A 


r 



L_ 


ABf|- (AB+Oos / '^y^) tan^\ 

A+X (A+Gos/\ / )tan^/\ 

^ c/4 c/4 


"1 



M-=r C ) 


(4.46) 
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where B = ^ 1 

On 

(a;:®4T=o ' 

Jj Jj 

M=0 


Oos^/N 


:/4 


variation of the low speed yawing moment 
coefficient due to roll rate at zero lift, 
given by 


^n 

(cf)c^=o 

]i=0 


A + 4 Cos 4\ 


tan 




3/4 


(4.47) 


where: 

X = distance from the G.G, to the aerodynamic 

centre, positive when the a.c. is aft of the 
C.G, and 0 is the m.A.C. 

A^n 

/ — , . R— ) rr effect of linear wing twist and is obtained 
6 

from Big, 25. 

0 = geometric twist in deg. 


AV 
( — £ ) 


a 


6 



effect of symmetric flap deflection (Ref, 5) 

streamwise flap deflection in deg. 

two dimensional lift effectiveness parameter 

(Ref, 5). 
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OSie tv/o methods mentioned above for the v^?ing contribU' 
tion could yield different values and the selection of the 
method to be used should probably depend on the type of 
airplane for which one is more applicable than the other. 


Using an elliptical lift distribution. Ref, 5 gives 
the following formula 


Jr 


8 


(4.48) 


(b) Vertical Tail ; She vertical tail contribution is estimated 
from the following relation (Ref, 5) 

- 2 (Vi) 0 

Py Py 

(ly Oosa+ Sin a) Cos a - 1^ Sina-j?;^) 


X Cy (4.49) 

Pv 


Ref, 1 gives an almost similar formula for the tail 
contribution 


G 

Py 



1 

Ma 

b 



0 0 -^ 


(4.50) 


^ • 4- YaW Rate Derivatives 

Yfhen a yaw rate, r, is super imposed on the airplane 
motion, the velocity field is significantly changed, Ror a 
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positive yaW rate, the velocity vector on the left v;ing 
increases and that on the right wing correspondingly decreases. 
This, in turn, increases the aerodynamic forces on the left 
wing and decrease them on the right wing. This lift distri- 
bution leads to a.n unsymmetrical trailing vortex sheet and 
sidewash at the tail. 


4.4.1 


) efficlent v/ith YaW Rate. C 


Y 


r 

4. 4. 1.1 This deriva/tive which normally is of minor importance 
can be calculated from the relation (Ref. 5) 


Oy -- «Y 

r r. 


2F«Y, 


V 


V 


- 2 


(l^ Cos a + Sin a) 
b 


G. 


(4.51) 


since the only significant contribution is from the vertical 
t ai 1 . 


4,4.2 Variation of Rolling Ifoment O o efficient with Yaw Rate, 

h 

r 

4.4.2. 1 This is another cross derivative, for a positive yaw 
rate, the asymmetric load distribution also produces a positive 
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rolling moment, which is proportional to original Cj^. Aspect 
ratio, sweep, taper ratio, dihedral and twist are all impor- 
tant. factors. The use of partial span hi^ lift devices also 
affects this derivative. The main contributions are from the 
wing and the vertical tail. If the horizontal tail is very- 
large, its contribution is significant and can be calculated 
in the same way as wing contribution and applying necessary 
corrections for effective sideslip. The vertical tail causes 
rolling moment proportional to its height above the roll axis, 

C. has little effect on the Dutch Roll damping but is 
r 

quite important to the spiral mode. Dor spiral stability, it 
is desirable that G-, "be as small a positive number as possible. 


4.4. 2. 2 Ref, 3 gives a very simple approximation for 0., as 

■^r 


1 “ 4 - 


(4.52) 


r 

Ror a better estimate, the contributions of various components 
are calculated as follows 





( 4 . 53 ) 


(a) Wing : The wing contribution is estimated from the 


following relation (Ref .5) 
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/A 0. 

(^) , 

M 


1 


(-^) JL 


Zag, 


1, 


r ^ (--^)e 


AG, 

“5 

6^ r 4 


(rad"’’) 


( 4 * 54 ) 


where: 


"1 

M 


variation of rolling moment coefficient due 
fco yaw rate at zero lift, given by 


^Oj=0~j 


f 1 I ) 


M 



G-, 


and (^) 


IfcO 


4B+2Gos/\, +2.. ■. 

.- c/4 tan /\ , \ 

( 


c/ 4 "^ ABf4GosA 


c/4 


8 


A + 2Gos /\^^^ 
■*■ A +iOos A 

fl 

X 

M =0 


tan A 


c/4 


j 


( 4 . 55 ) 


, G,-0 Variation of low speed rolling 


mom ent 


coefficient due to yaw rate at zero lift 
and is obtained from fig . 26 as a function 
of aspect ratio, sweep of the quarter 
chord and taper ratio. 
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increment in 0., due to dihedral, given hy 
■^r 


^ Tt’ASinA^^^ 

T2 "1. + 4 CosA 


(4.56) 


-) = increment in due to wing twist and is 
obtained from J'ig«27 


\ 4 


effect of symmetric flap deflection (Ref. 5) 


and a„ and 6' are defined earlier, 

U ^ X 

Refs.l*, and 4 give the value of (p ; — •) as a function 

^1 

of aspect ratio, sweep and taper ratio in the form of charts. 

(b) Vertical Tail ; The vertical tail contribution can be 
estimated from the relation (Ref ,5) 


- 2 (|) (|■).CY 


(l.r Cos a + Sin a) (z-rr Cos a ~ 1.,. Sin a) 


(4.57) 



(c) Ho r i zo nt al T ail : Ref, 2 gives a relation for estimating 


the horizontal tail contribution 


G. 


1 , 


-H 






Ph 


JL 


bS 


(4.58) 


Q 

v/here ( l^,)^ is calculated as for the wing from (4,54) and 
(4.55). 


4.4.3 Variation of Yawing Moment Coefficient with Yav>^ Rate , 

r 

4. 4 . 3.1 Ihis is also called the 'damping-in-yaV'? ' derivative, 

since it represents the resistance of the airplane to yawing. 

Ihe unsyrametric aerodynamic forces cause an increase in the 

profile and induced drag on the left wing with a corresponding 

decrease on the ri^-t wing, .This produces a negative yawing 

moment. The effect depends on aspect ratio, taper ratio and 

sweep back. The side force on the vertical tail also produces 

a negative yawing moment, .The main contributions are from the 

wing and vertical tail, C is the main contributor to the 

damping of the Dutch Roll mode and also plays a significant 

role in determining spiral stability. It is a very important 

parameter in airplane lateral stability characteristics, Ror 

best effects in each of these modes, large negative values of 

C are desirable. 
n„ 



4. 4. 3. 2 The contributions of wing and tail are estimated as 
follows 




+ G 



(4.59) 


(a) Win^: The wing contribution is estimated from the relation 



(4.60) 


where: 

G 





a factor accounting for the variation of induced 
drag and is obtained from Pig. 28. 

a factor accounting for the variation of profile 
drag and is obtained from Pig, 29. 




wing profile drag coefficient given by 

(4.61) 


(^) Vertical Tail : The vertical tail contribution is estimated 
from the relation (Ref .5) 



2 (V G, 


(l^ Gos a + Sin a)‘ 


G. 


■ b' 


Y. 


2 


(4.62) 



MD DISCUSSION 


5.1 The methods of estimation of lateral stability deriva- 
tives have been applied to a sample case. The Gessna-182 
airplane has been considered for this purpose. The mass 
characteristics of the airplane are given in Table 1. The 
geometric characteristics and other allied constants which 
are essential for the estimation of stability derivatives 
are given in Tables 2 and 3. These constants do not vary 
with the flight conditions. Table 2 gives the variation of 
fuselage cross-section area along the body length. 

The two flight conditions considered, approach at 
73 Diph and cruise at 150 mph at 5000 ft, are given in Table 4. 
Table 5 gives the list of parameters dependent on flight 
conditions. 

Based on the methods of Section 4, lateral stability 
derivatives have been estimated for the two flight conditions. 
The values of derivatives along ¥/ith the contribution of ma;ior 
components as estimated by the methods of Refs.1 to 5 are 
presented in Table 6(a) and (b). The values estimated from 
Ref. 5 ar-e complete and compare very favourably with the values 
provided by the manufacturers which are also presented in 
the same tables. Table 7 mahes a comparative study of the 
estimated derivatives and their values supplied by the 
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manuf acturersj in tlia'fcj it gives tlie percentage error from 
these values. 

5.2 In some respects, the methods of estimation given in 
Refs. 1 to 4 are not complete, in that, all the contrihutions 
are not accounted for, a discussion of these methods and 
values obtained from than follows, 

5.2.1 Cy ; Refs. 1 to 4 do not provide any theoretical 

relations for estimation of this derivative. All of than 
suggest use of wind tunnel test data for the design in question. 
According to Ref .4, interference effects are so large that a 
generalized formula would not be completely satisfactory. The 
method of Ref, 5, a very recent publication, Includes inter- 
ference effects based on experimental results and hence probably 
gives the most accurate theoretical estimate. 

5.2.2 G. : The methods of Refs, 1, 2 and 4 give a slightly 

P 

higher negative value of wing contribution than that of Ref .5, 

This is because of the variation in the values of the factors. 

Only Ref .5 gives a method for estimation of horizontal tail 

contribution. The vertical tail contribution is almost the 

same in all cases, Tiie overall value of G-, from the methods 

P 

of Refs. 1 to 4 is slightly higher (but compares more favour- 
ably with that of the manufactures) and from that of Ref. 5 is 
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5.2«3 0 5 Ref, 1, 2 and 4 have given relations for estima- 

ting wing contribution. This contribution being proportional 
to square of lift coefficient, is considerable for the approach 
configuration. If this contribution is considered, the value 
of G becomes too lange compared to manufacturer’s value. 

However, it is negligible for the cruise condition. The other 
references neglect this contribution entirely. The body contri- 
bution fram Refs, 1, 2 and 4 is a much higher negative value 
than that from Ref. 5. The method of Ref .3 gives a very low nega- 
tive value for this contribution. The tail contribution is 
almost same by all the methods. The overall value of 0 

n 

p 

estimated from Refs, 1 and 4 appears to agree quite v^ell with 
manufacturer’s value for approach configuration only. Ror the 
cruise configur-ation, these two values ai'e too low. The value 
from Ref. 2 is too lov\f and from Ref, 3 is too hi^ for both 
configurations. The value given by Ref, 5 is more consistent 
for both the configurations, 

5.2.4 Gy 5 As mentioned before, this is not a very important 

derivative. Refs. 2 and 3 neglect this derivative. The methods 
of Refs. 1, 4 and 5 give small positive values for this deri- 
vative since the angle of attack is more. However, the mai^n- 
facturers give a negative value for this derivative. The methods 
of Refs. 1 and 4 take the effects of sidewash into account. 
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Ref, 6 quoting USiR 'DAfCOM*;, gives a method which gives a 

value of Oy much in agreement v>/ith manufacturer's value. 

P 

But since this derivative is not important, the value given 
by Ref, 4 only is considered in the final equations. 


5.2.5 : Refs, 1, 2 and 5 give the same formula and charts 

for the wing contribution. The value from Ref ,5 also is the 
same. The method of Ref. 4 also gives almost the same value 
for this contribution. Refs, 1 to 5 do not give any relation 
for the horizontal tail or vertical tail contributions. The 
methods of Refs. 4 and 5 give almost the same value for these 
contributions except that the relation for vertical tall con- 
tribution given in Ref, 5 is sli,^tly modified in Ref. 4, The 

overall value of Gn compares well v^?ith the manufacturer's 

P 

value. 


5 . 2.6 




The methods of Refs, 1, 2 and 4 give a smaller 


negative value than that of Ref. 3 or Ref. 5. The tail contri- 
bution is almost the same by all the methods. Correspondingly, 

the overall value of 0 by methods of Refs.1, 2 and 4 is 

^P 

less negative and that of Refs. 3 and 5 more negative than 
the manufacturer's value. Ref, 6 gives a formula for 0^ from 

from WAGA~TR 258? which gives an 



<il7i 


NAGA-TN 1581 and for C 
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overall value which a.grees very favourably with manufacturer 's 
value. However, value given by Ref, 5 has been used in the 
solution of equations of motion. 


5,2,7 0^ : Ihe value of 0^ by the method of Ref.1, neglec- 

ting the effects of sidev'/ash compare very favourably with 
manufacturer's values. However, the value by Ref, 5 is consi- 
dered for the solution of equations of motion. Refs. 2 and 
5 neglect this derivative entirely. 


5.2,8 0-, • Ibe value of C by all the methods except Ref ,3 

are almost same, Ror the approach configuration this value 
is as much as 307, more than the manufacturer's value. Rof 
the cruise configuration, the agreement is much better, The 
value by Ref ,3 is too rough an estimate. 


5,2,9 i All the references give the same relation for 

this derivative, The estimated value compares favourably with 
the manufacturer's value, 

5,3 A consolidated list of estimated derivatives and those 
supplied by the manufacturers lia,s been given in Table 8, 

These two sets of va^lues for each of the tYW flight configu- 
rations have been used in the solution of equations of motion 
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for an Initial disturbance of 10° in sideslip angle. Typical 
input/output for the digital computer programmes of Appendices 
’0-' and 'E' are shown in Appendix 'E*, The resulting response 
of the aiiplane in terms of roll, yaw and sideslip angles 
and roll and yaw ra,tes are presented in Appendix 'O’, 

5.4 Computer programmes were developed for solution of the 
' equations of motion, both by Laplace transform method 
(Appendices ’A’ and ’G’^^Eourth order Runge-Kutta method 
(Appendices 'D* end *E*)* Ihe former method computes the roots 
of the characteristic equation in addition to the time history 
of airplane motion, v\?hereas, the latter computes only the time 
history, but has the advantage of being simpler and quicker. 
The results obtained by the two methods agree very closely 
for the time increments used (.05 sec) in the numerical integ- 
ration, The roots of characteristic equation obtained from 
the Laplace transform method enable the estimation of the 
period and damping characteristics. These characteristics for 
the two fli^t conditions are given in Table 9(a) and (b). 



6> comusioiT 


6.1 Kie methods of estimation of lateral stahility deriva- 
tives from Ref* 5 appear to be more appropriate, since, they 
are complete and the derivatives estimated from these methods 
for the sample case of 0essna-182B compare very well with 
those provided by the manufacturers for both approaph and 
cruise configurations. ' 

6.2 The application of the methods of this report to the 

above sample case predicts dynamically stable lateral airplane 
motions. In the spiral mode, time to damp to half 

amplitude, is 125 secs, and 62 secs, for the two flight condi- 
tions as against .53.5 secs, and 38,5 secs, predicted by the 
manufacturers values. T his shows a tendency towards spiral 
unstability which seems to agree with the actual case. 

6.3 The va 2 ;iation in and damping from the manufacturers 

values for the other two modes, namely, Rolling and hutch 
Roll, is within 5%' and 10% respectively. 



THIS PA.GE IS LE]?T BLiHK IHTENTIONAILY 
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'IHIS PA.G-E IS lEIT BLMK lEDENTIONAlII. 
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SPECIAL FOTi^TIOHS USED IN SECTION 


1 


If two roots A^. ^2 Of the oharaot eristic equation 

are comply, then we can write 


A,, A, 


R + i I 


(B.1) 


Purther, we can write the power of as follows 


- R k 

/\ = (R + X I) 


1 


_ pk ^ , „k~l k,(k-1) k“2 2 

-R+kR (11)+--^ i_R (ii) ....+ (1 


(il) 


= R*" - k i2 + jEA3EzlL jR-2)-(krXi. jj(k-4) ^4 


I^- 


k.(k-1) .(k-2) R-*3 3 

+ i(RR I .R, . . ) 


" \ i ^k 


(B.2) 


We can now substitute for in the expression for 

the coefficient A. 


= £ a? (R. . + il .) 

‘ j;=0 7-D 5-.1 ' 
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+ Ur 


(6AR^ + 5BR^ + 4CR^ + 3PR^ + 2ER^ ) 

+ i(6AI^ + 5Bi^ 401^ ^ ^ 2EI^) 
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which can he written as 

. = ^ _ .^1 ^2 + y ; 

1 + i “ V 5 - ' — 


A. 


^2 + y? 


^2 + ^2 


- R, + i I, 
A ^A 


(B.4) 


A 2 will be the complex conjugate of A^ and will be 

"^2 = \ i (B.5) 

Similarly 


®2 - t. i Ib 


^1’ ^2 “ i ^ 


(B.6) 

(B.7) 


If the other pair of roots is also complex, then 
exactly as above 


t= 

"1 

+ I 

I ' 

(B, 8 ) 

B.., B. 

3’ 4 = 


+ i 

I' 

Xb 

(B.9) 

G 3 , = 

I'd 

+ i 

I* 

G 

(ii.10) 
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I*^-? PI^0-^AM.M£ evaluates THE COcFF I G I ENTS OF CHAKACTcR I S T IC 

MOTION OF AN AIHPLANF USING MASS, GEOMETRIC 
AsKUiJY NAHiG PARAH-T=RS, SOLVES THE QUARTIC cQUATON, SOLVES 
^0UATICN:3 of MOTION BY LAPLACE TRANSFORM METHOD TO PREDICT 
fIM- HISiORY OF ROLL, YAW AND SIDESLIP ANGLES AND ROLL AND 
RA TSS . 
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CL 3 =. 5 -CL=s=TAN{ GAMMA) 

•;\K=( (XKO^COSI ETA) )**2 + ( ZK0«SIN(ETA) ) ) / { SP AN**2 ) 

ZZK = ( ( ZKO-CUS ( tTA ) ) + ( XKG^SI N { ET A ) ) ) / ( SP AN«*2 ) 

CZK = ( ( 2KO-X*2-aKO= 5‘*2 ) *S IrK ETA) *COS( ETA ) ) / { SPAN^*2) 
-ki=xzk/x.-;k 

AK2=^XZK/Z2K 

TaU=AMA SS/(RHO=i=AR£A=}‘VeL) 

F'L R I =2 • I ^TOU 
T HN C 9 '*TOU 
AM!J = AMA SS / ( RH j«AR£ A*SP AN ) 

SL0=.5*AMU/XXK*CLB 
SLP=0.2 5*CLP/XXK 
SLR= .25*CLR/XXK 
SlC- .5* AMu/XXK^CLC 
5NB=:.5*AM!j/ZZK*CNR 
SNP = 0.2£=5‘GNP/ZZK 
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S;^R=G. 2 5*CNR/ZZK 
SGC=.5«4HU/ZZK^CNC 
Y B = . 5 ^ CY ;;■ 

YP=.25*CYP/AMU 
YR= .25*CYR/ami) 

YL=.5«CYC 

^ - SL P— b :iR + aK i’i*SNP+ AK 2*SLR. 

‘^^=SLP*SMri-SLR*SNP 
P- =SLB* SNH-SLR*SNF. 

P 1 , = S L P * S f'J S - S L B * S ! J P 
p:- = AKi * SfJb-SLB 


Pc = «K2*SL!:'-Si\IB 

P7=-PZ=^ Vb + P3*YP + P4*YR-PA 

Aj?icNSiOi Y(t,) 

( ~ A K i A K 1 

^ ( 3 ) = P:,-a { 3 )=f'-YB 
'^‘5=~^I*Y0+P2+P5*YP+?6*YR-P6 

1 ( -1 ) = p 5 «c L + P6 LG + P 7 

A { 6 ) =PB =». 3^CL + P4*C! G 
PPviGT 3C. ; 


PBIBT 
'^RlhT 
PRI ■'!! 
PRIHT 

PR I . at 


- CO;. , WLI T , X I KO , Z IKO , SPA:j , ARc A 
:• CO 3 , V r L M » RHO , C L , 4 L F AD , G AMM A D 


? »CYB»CYP,CYR »CYC,CLB,ClP,CLR,CLC,Ci\t3,Cf4P}Cf-jR,CNC 
.',PHln,OPHIU,PSlO.DPSIO.BTun i>n-fUWK,UWL 


r f ) 


PRINT PC2- 
•PRINT 3, i 
PRINT 5^-_2 
I i4T . Co _ 
PRINT 3027 


PR T NT -L , { A ( I ) , I =2 , 6 ) 

.'~u=PHiD=«'A(2) 

A^=PHIG>>A{3)+DPHIQ*A(Z ) 

;^£; = KHia=«‘A { 4 )-BTAa=^P5 + DPHla* (-A ( 2 ) *YB+ AK2*SLR-SNR ) 

- up b I ili* {AKl*Sr<R-SLR) +SLC-Si4C*AKi 

cTd PSIO*P3*CLG-BTAO=!'PO+OPHIO>}^(P6«YR-P6-AK2^ 

'^-^^’^^0*{-P&«YR+P5+AKl*SNR*YB~SLR^Ye)-SLC*{SNR + YB ) + 
_bNo4£(AK i*YB + SLR )-YC*P5+OPHlO=!'SNR’i=YB 
Ax.= ^ PHI c=!«p^-PSIO*P3 + DPHiO*P6-DPb]a*PB )=i^CLG + SLC«( -^NB-SNB* 

: YR+SNR«VBI+SNC=«=(SLB«YR-SL8-SLR«YB)-YC=f'P3 
A5 = CLG* (-0LC«Sr4B + SNC*SLB ) 

PO=PSIO>!'A{2) 

B> = PSiiJ^A(3)+DPSIO*A{2 ) 

B2 = ^S!0=»A{4)-3TA0’i'P6-DPHl0«(AK2=i=SLP-SiMP)4DPSI0=!‘(-A(2)=J=YB + AKi=!< 

.SiiP-SLP )-SLC*AK2 + SNC + 

u3=-.5*FHIO*P6«CL + PSIO^{ . 5*P5*CL + P7 I-BTAO^PA'+OPHI 0* {-P6*YF + 
iAK2«SLP ^i'YD-SfJP^YB ) +DPS 10* { P5*YP-AK1 *SNP*Y8 + SLP*YB ) + 
2bLC*{AK2*YB+SNP)-SNC*(SLP+YB)- YC*P6 

PSI0*P3-DPHI0*P6+DPS10*P5)*CL + SLC*(SNB*YP-SNP*YB 
* + SNC*{SLP*YB-SLB*YP)-YC*P4-,. .5*PHI0*P4*CL 

B5=,5*{ SLC*SNB-SNC*SL3 )*CL 

C u ^ b T A I J ^ A { 2 ) 

C.= .5*PFi J*A(2)*CL + PSIO*A(2)*CLG+GTAO*P* + UPHiO*A(2)*YP-DP<?Tn* 

iArZ )*( YR-_,, )+YC*A(2) urbiu^ 

Ct;=:.5'i'PFIuvpj.T:CL+PSIO*Pl*CLG+6TA'j*P2+DPHIO*(.5=i'CL*A(2)-AK2*SLP* 



nl 


i 


■K)5, 


.000 


3i00 




J^Ou 


4YR+Ai<2*SLP 

Of* c cho-D^ 2*P2*CLG+DPHIO=<‘(-AK£*SLP*CLG+SNP=i'CLG+AK2>f-SLR* 

“^,o!Yf:!40rr:Or"^'^'^*^^''*^’^^'^"^'^2=«'CLG) + Sr4C«(-SLP*YR+SLP4- 
-J 5L \*Yh AK *CL*. 5+CLG ) + YC=!‘P2+u . 5=i'CL*PH I O’i'PZ 
C4=SLC* ';P=i'CLG-5NR*.5>s<CL) + SNC*(SLR«CL«.5-SLP*CLG) 

CALL QUARK A, Y) 


'"K i I = i , f 

T(!ALF{ i )=THN/Y( I ) 

IF {Y{i). Q.j.<i) on TO 2 
F^R< n=Per.T/ABS( Y{ 5) ) 

I F { Y { f) .ErO.O. 0 ) GO TO H 
fKR { C) = F = RT/A8S{ Y( 7) ) 
COUTIOUE 


PKloT I (PcRU ),I=1»2) , (TH4LF( I ),I = i,4) 

FjRMiT ( /5.K6{ F20.4) ) 
aA ( 5 ) = a 5/ .. { 6 ) 

AA{o}=( A4-A5«A( 5)/A(6) )/A(A) 

R3(3)=e5/', {5) 


Ce.(6) = ( B4-B5>^A(5 )/A(6) )/A{6) 
CC( 5)=C A/-. (6) 


AAY { 5) = AAS.5)/T0U 
U8Y ( 3 ) = 63? 5 ) /TOO 
IF{Y(5) .fK.U.U) GO TO i5 
IF( Y(7) .0 :. j. .^) GO TO i:^5 
F* vi s j I “ A f 4 

AU( I ) = A£+Y( I )* {A4 + Y{ I ) =»= { A3+ Y ( I ) * ( A2 + Y{ I ) =i= { Ai + Y ( I )*A0) ) ) ) 
-;(I)-Y(n=f={2.*A{6)+Y(I)4=(j.«:A(5) + Y(I)^{4,=f*A{4)+Y(I)^(5.*A(3) + 
_ Y ( I ) 4^6 . =("4 ( 2 ) ) ) ) ) 

AA( i ) = iC{ I )/E( I ) 

BlM n=R5+Y { I )* (B4 + Y( I ) *(B3 + Y( I )«=( o2 + Y{ I ) * ( Bi +Y ( I ) ^Bu ) ) ) ) 

BB{ i )=PC{ ; )/E { I ) 

CD( I )=Y (I )>!'(C4+Y( I )*(C3 + Y{ I ) * ( C2 + Y ( I ) * ( C 1 + Y { I ) *C0 ) ) ) ) 

CC{ n=CC( I )/£(] ) 

A AY ( I )=AA ( I )❖¥(! )/TOU 
BfcY( I ) = eB(I )>!'¥( I )/T0U 
COUTIU5JE 
PRIUT 3',’ j 


FOR -KTJ //-0Xf=^'AA{i IlXf*AA{2) = '^jiiXt*AA(3)=^,'» iX»’i‘AA(4 ) = *, 

i„X,*AA (5) = *, 11X,*AA{6)==^) 

PRI4T BUBA, ( AA{ I) , 1=1 ,6) 

FOR 4 AT { /4X »6 ( 4X ) FS . 5 } 4X ) ) 

PRIUT 3100 


FORMAT ( /'/. 0X,=«=BB{i.)=*,llX,*BB(2)=^^,iiX,*BB(3)=*,liX,^B3{4)=#, 
i. ilX,*3e( :>)=*, llXf*BB( 6)=«) 

PR I UT 2’! 3: = f ( 3B { I ) j I =1 }6 ) 


PftlOT 32 uO 

format { // ‘fKt^CC{i)=^, I 

i^X ,^CC (5 )=* ) 


X,^CC(2)-*,iiX,^CC(3)=*fiiXt=i'CC{4) = *, 


PRIMT 2. : yKCC{ n,I = it5) 

C'RINT 3000 


FORMAT ( / / ,, jX , *AA Y{ 1 ) = =5^ , 1 1 X , =!'AAY{ 


= *j ilX,=j'AAY ( 3 ) =*, iiXf 



7S' 


dO 


9:' 

lOiS 


1*aAY{ 4) =*jliX,=«‘AAY(5) = >J=) 

P;^INT 30PA,{AAY(1),I = 1,5) 

PRINT i 

F.jRMaT I // - dX,*BBY{ 1 )=«,iiX,«88Y(2 ) =*, iiX , *BBY { 3 ) = =^'t 
)==^=, Il,X,=«=8BY (5)=*) 

PRINT 30i:., (BeY( 1 ) ,1 = 1,5) 

PRINT 3C27 
PRINT 3C. 

PRINT 3CZ;; 

PRINT 3£:,. 

PRINT 3;>::;7 

PRINT 6C,(AA{ 1),Y{I),I=1,4),AA(5),AA(6) 

Fur MAT ( // ;X . = *,4{ F10.6 , «tXP{ SIGMA*, FlO. 6, *) + *)/ 

i.5X, Fr„ . £, i=SlGMA+^, FI.' , 6) 

PRINT ej, {BB(I),Y(I),I=i,4),oB{5),BB(6) 

PRI JT SO,(CC( I ) ,Y( I ), 1=1,4) ,CC{5) 

PRIRS 90, (AAYII ) ,Y( I ) ,1=1,4) , AAY( 5) 

PRINT 9 ,{BRY(I),Y{I) ,I = i,4),RBY(5) 

format ( //bX,* = *,4I FI U 5, «i;XP( SIGMA*, Fi - .6,* ) + ^ ) / 

^5X, Fj 0. 6) 

T = 0 , 0 

00 LIK:=^,.2>.= 

UiMrOSI C . AB( I- ) ,BA (1C),BC( iu ) , bAAY ( i. ;. ) , ABB Y I'i J ) , ZD ( i> ) 
SIG'-’A = T/Tf^U 
Sum =0. 

SUM P = -, ... 

s:-ip =-. . c 

5: '.UP =0,0 
TOT =0.0 

00 i.i. l=i,4 
p = S1GmA ’i'YI I ) 

ZD( I ) = OXP(P) 

641 I )=AA( i j^i^ZUI I ) 

B^AYI I ) =A^Y( I )*ZD( I ) 

£B( I )=3e{ : )*ZD( I ) 

A30Yi 1 ) =8L-Y{ I )*Z0( I ) 

3C{ I )=CC( i )*ZD{ I ) 

5UM = SUH-+BA( I ) 

5UMP = SUf'P + SA.AY{ I ) 

SHP^SMt’+AGII ) 

SMUP=SMLP+A38Y( I ) 

T0T = TuT4B:,( I ) 

CONTINUE 

PHI =5UM+AA ( 5 )*SIGMA+AA (6) 

PSI=SMP40&{5)*SIGMA+8B (6) 
r3£T.a = TUT+CC{ 5) 

R0LVEL = SU 'P + AAYI 5) 
yMWVrL=Srv10P + BBY(5) 

PHIDILI Kl ) = PHI*13:). /pi 
PSI CMLI K£ )=PSI*18d ./PI 
6tTi.0(L IK.; )=D£TA=5‘iaO./PI 
ROVELDI LIKE )=RQLV£L*ISC./PI 
YaVELDI LIKE )=YAWV tL*l80./PI 
m LIKE )=T 
T=T+0.'j5 



7*7 
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Sv CONTINUE 
GU TO 65 

' Rl=Y{il 

= Y ( I ) ( 5 ) 

R3=y (1) 5-3 .*Y( i )«Y( 5 )*❖; 

SlII M!! 5) )i*2 + Y(5)>i‘*4 

AI ~ = Y{ 5 1 )=i= 4 : 3 «Y ( 5 )^^Z + 5 ,^Y(l ) *Y( 5 ) #*4 

1 1 2 = Z . V { 5 ) ^Y ( i ) 

a1 3 = 5 .=S« Y ( ^ ) =s<*2«Y { 5 )-Y ( 5 ) 454:3 

AI+ = *4..4(Y(1)4434Y(5)-Y(1)4Y(5)443) 

''■''J-=*'R6-+A245Ri + A34R2+A44Rl+A5 

AY 1 =A 04 AI 5 + Ai 4 Al 4 +A 24 A I 3 + A 3454 I 2 + As- 45 Ar i 

YvYOtsui’^i-Ij’tN^!* Ri*5.»il5>« S2+2.*A(6>» 

KA = { AX J,*.<v + A Yi *y 2 ) / 2 
aIA = ( 2 2 * 'iY. -A'< 145 Y 2 ) /z 
i3x;i=:aQ#R5^.g|4,r^^^g2 4R3+e34H2 + 8445R ' +B5 
NYl = 304 ^ 4 1 ;. + 61441 6 ’ + 324 <aI 3 +B 34 ^A 1 2 + 64*411 
60 =(BXi4.'2 + BYi*Y2)/Z 
316 = { / 2 * ,Y.-3Xi*Y2 ) /Z 
Cx.=Cu*R5+Ci*R4+C2*R3+C3*R2+C4*Rl 
Y - =C J*4i : + C i* A 1 4+C2 *A 1 3 + C3*A 1 2+C4*A 1 i 
RC =(C< i 4.,2 + CYi*Y 2 )/Z 

CIC ={ A 2 *CYA.-CXi*Y 2 ) /Z 
aka = 2 . * SORT ( RA** 2 +AI A** 2 ) 

3X3 = 2 . 4 SQRj ( RR442 + B j ^2 ) 

CKC = a .*SQRT ( RC** 2 +CIC**'’ ) 

'^'\CT^N{RA »A 1 -\ jWA ) 

CALL ARCTAN(RB, BIQ, WB) 

CALL aRCTaNCRC ,CIC , WC ) 

CAL L aR CT {Y(i),y{ 5 ),WAA) 

YAKA = 5 CRT{ Y( i ) ** 2 +Y{ 5 ) **.? ) /TOU 
AKAA=AK A 4 I'AK A 
AKA ri = y K B* YAKA 

IF( Y{ 7 ) .N , 0.0 ) GQ TU 125 
00 5 : I =;, , 4 

AD( n=A 5 + Y{I)*{A 4 +Y(I) =S=(A 3 + Y(n 4 (A 2 + Y{I )*(A 1 + Y(I ) 4 Aa)) )1 
AA( i ) = AC{ I )/ 2 ( I ) 

BDU ) =B 5+Y ( I ) 4 ( B4+ Y ( I ) 4 ( B3 + Y ( I ) 4 ( b24Y { I ) 4 { B1 + Y ( I ) *B0 ) ) ) ! 
;-^;;;=jfi)^jc^+Y{i)4(c3+Y(i)4{C2+y{i)4(Ci+Y(iUco)))) 

CCd )=CC( x)/z( I ) 

AAY( I )=AA ! I )*Y{ I )/TOU 
G 6 Y{I) = eE(I)*Y(n/T 0 U 
CuNTINUg 
PRINT 55 > 

^ 7X,4AA( 1 ) t AA (2 )=4, 17X|*AA{ 3 ) = *, 1 iX , *AA ( 4 ) =* , I a X , 

R i\ I i 4 1 A 3 ■ u f K A f A I A } ( A { I ) f I = 3 I 6 ) 

FORXlAT{ /l 'X,F3.5»*+/-I*,F3.5t9X»4(5XjF3.5j5X) ) 


Ri 



80 


i -'i > 36 Ju 

. 7X , «hq ( • ) / -- « j, 

37 L 0 ■^'^’^"'"'’'■'''■'> 1 . 1 = 3 , 6) 

{//^7X,«cC(i ),rc ^., 

' ' ’-^■'>*BC{3)=»,i;x,*cc(4i.»,ux, 
, >'C,cic,(cc(n,i,i,5, 

’ . 6 .'' f * Y ( 3 ) = * .V .4 

. 6] (/ \f I ' 

AKT- 3p:;3’“ ) 

^ I -> i ' ! ; V ■, <« i { / / 1 X , R p Y / ^ I j, . 

J ' 66 f (d3Y{ I ) , i = J’r 1^^»*B3Y(5)=7 |c) 

J 30;.? ^ ^ 

ORI 7? 3 i 

'KI :T -. 2-~ 

print 3 : 

PRI^VT 23^7 

PRINT i 9 . r AK^ , Y { • ) . Y f c , .,, ,,. 

iM. .ras«r(//:x,».....;.;:i;;'?''‘'J'”;j;i),i=i,4),4A(5),A4(6) 

FORMAT ( //: , , . , . c 

2F10.6) 

_ 4*^.0 FuI^M^T ( //fX}^, . - 

. wnt f -''«r.vO),y,5;,':i;SaI!I,:a^ ; . 

i045 T = 0.0 -’ '■'’’^'■’'“''-"«''»BX(n, Yd), 1 = 3,4), BUYIs! 

i^O i,0^-j L i C K = i j ^ 3 u 
SI-5M4=T/T ;iJ 
SUM =.;. 

SUHP=D, C 
SMP . V 
SMUP =). 

TOT= ^ . .- 

00 106 C 1=3,4 

P=S I GMA *Y ( I) 

ZO{l|=tXPlP) 

3a( I )=A A( I )^ZD( 1) 

eA4 Y( I ) =a; Y { i ) *Z0( I ) 

At ( 1 ) = *:■ Ell )*ZU ( I ) 

ABBY ( I } =3 ,Y{ 1 ) ^12 ( I } 

BC ( i ) =Y C I * ) ^ Z u ( I ) 

sy,M=su^* +e . ( I ) 

SUWP = SiJFP+.Ba4,Y( I ) 

S«P = SmP 4a -; ( I ) 

SMljp = SMLP + A3BY{ I ) 

TOT = TOT +63 ( I ) 


SI 


i05;. Cu^MTlNUE 

CUSA=C0 S{ iIGMA»Y{ 5 ) + WA ) 

COSB=CnS( SIGfiA*Y(«>)+Wd ) 
u U S C — C '■ j S ( .i I w A ^ Y ( D ) + W C ) 

CUS?=CO S { SIG'^A^Y ( 5 ) + WA+WAA ) 

C GS “ C tJ S ( o I G '‘1 A ^ Y { Z ‘ ) + b + W A A ) 

GXPR = cX F( jIGMa*Y( 1 ) ) 

PHI =AKA=^‘c..PR*CrjSA + SUH+AA( 5 )*SIGMa + Aa( 6 ) 

PS! =BK3=5't vPR^CUSB + Smp + bb ( 5 ) *S I GMA+Bb ( 6 ) 
e;TA=CKC* iXPR^COSC + TGT+CC ( 5 ) 

RULypL= AK-4A*'iXPR=«^CGSP+SUMp + AAY( 5) 

YAW V = L= AK B^S^-XPR^CGSR + SMUP + BbYC 5 ) 

PriID(Ll CK J^PHI^ieo./PI 
u{LI CK) = PSI*I30./P1 
B. TA0{L ICK) =B.:TA*iR ./PI 
ROySLOt LIuK )=R0LyEL*j.30. /PI 
Y„ycLD{ LI gK) =YAW ypL*l30./Pl 
. Tj {LICK, )= ; 

T=T+ . 5 
„;-=4 Cj-ni;iUE 
Gj to e£ 

1*5 i)U :,0b i = i,L 

Ag( i ) = A5+/( 1 )-{AA-+Y{ I ) *( A3+Y( I ) ^ ( A2 + Y ( I ) * ( A*4 Y ( I )«A0) ) ) ) 

. { ! )=Y ( I )l'( 2.*A( 6)+Y( 1 )*{ 3.#A( 5) + Y{ I )’i'{4.*A (R )+Y( 1 )*{ 3. *A (3) + 

_ Y ( ! ) *6 . =^A { 2 ) ) ) ) ) 

'\,4 I )=AC( n / ; ( I ) 

BP( 1 )=B5+Y( 1)*(B A+Y( I ) ❖(e3+Y( I ) * { B 2+Y ( I ) « ( +Y { I) *BU ) ) ) ) 

3P( I )=^^C{ I )/={ I) 

CL)( I )=Y {1 1 )*(C3+Y( n«(C2+Y( I ) * (Cl + Y (1 ) ^S'CO ) ) ) > 

CC( 1 )=CC(I )/c( I ) 

AAY? I )=5Ain*Y(I )/TnU 
6bY( I )=e0 { I )*Y( I ) /TOU 
.05 Cj'IIlHUS 
.v2? R.=Y(;) 

R.=Y(3M=!',-Y(7)**2 

Rr = Y{3 ) -3.*Y( 2 )*Y{7)*^2 

R4=Y{3 ) =t*--6.*{ Y(3)=)=Y( 7) ) **2 + Y ( 7 ) 

R5 = Y ( 5 ) '-1.V . =f=Y { 3 ) ^^3 *Y { 7 ) **2+5 .*Y ( 3 ) *Y { 7 ) **4 

Ail=Y{7 ) 

Al2=2.*V(n*Y(3) 

AI3=3.*Y{ . )**2*Y{ 7 )-Y{ 7) **3 

{Y(3)**3*Y(7)-Y(3)*Y{7)**3) 

Ai5={5. *Y(3)**4-1 '.*(Y (3 )*Y{7 ) ) **2 + Y { 7 ) **4 ) *Y { 7 ) 

AX I =AO*R5+Ai*R^+A2*Ri.+ A '* ^2+A4*Ri. +A5 
AYi=Au* AI i+A I*AI4+A2*AI3+ A3*AI 2 + A4*A 1 1 

;U=b.*A{2)* R2+5.*A{3)* R4+4.*A(4)* R3+3.*A(5)* R2+2.*A(6)* Ri 

yI= 6.*A {i )*AI5 + 5 .*A(3 ) *AI4+4.*A(4)*AI3 + 3.*A{ 5)*AI2+2.*A{6 )*AI 1 
Z = X i_**.3 +Y. **2 
R A D = ( A ■< 1 * •. 2 + A Y i. * Y ; ■ ) / Z 
41AD= ( X 2* ';Y*-AXi*Y2) /Z 

eXl=P0*B5+Bl*R4+32*R3+B3*R2+B4*R*+B5 
<;Yl = 3u*AI 3+Bl*AI A-+B2*A IX+33*AI2 + uA.*AI 1 
RbD= ( 3>; 1* -2 + 8 Y l*Y2 ) / Z 
8I30={X2*iYl-BXi*Y2) /Z 
CXi =C3*B5 +Ci*R4+C2*R3+C3*R2+C4*Ri 
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■15 


:90 


. 0-55 


i-Yl 54 .ri * A T 

nCD = lcxi*x2+cyi*y2i/7 
CiCL)=(XJ*CY.-cxi»Y 2 ,/z 

r'r ■ o- = 1 CO** ’ 1 

AIAD, WAD) 

r^M f^^l''^^«B 0 ,BIBD,WBD) 

.^Kpu-5CR 1 ( Y(i )*«2 + Y{7 )**? ) /TQU 

'^"V'^CT.ri(Y{,),Y{7),wAA5; 

auA 4 J=m K Al;*YmKA0 
AKABn = BK6>.'*YAK4D 

s^ivlMT ^x 5 - ,iUX,i(Fd.5,5X) ) 

F.3R 1AT{ // _7X,*RB( ■ ) ,3R 7 . * 

iii’'E5:Y?r;-“----,3,=Y, 

r jRvAM 5..<,2{5X,Fe.5,5A)} 

PKliu* 2 u '7 

PHI XT 5C , 

PHliT 1 C 2 : 


PRI.mT 
PR I NT 


300 i. 

2 u;_? 


PR I JT 
PRI-,'T 


I tAKA,Y( ij , y{5),wA,AKAD, 
RI-Ji^9C, CKC,Y(1 ) ,Y{5 ),WC,CKCD, 


01^ 


J 'xT 3 C 0 ^ 


PRINT 


AKAA,Y(i) , Y(5) ,WA,WAA, 

J.. H M ^ ^ 2 ( 3X, F I 

^ T fZ L4 A. C -- r - ? — \ t r ^ 

iH+»Pl A | 2 H)Hf*/) 

fZC^XfFj, 

f 2 H ) +/ ) 


Y ^} tY( n, W6U,B8(3) ,BB(o) 
»Y(7) ,WCD,CC(5) 
AKAAD,Y(3 ) ,y(7 ) 


' i f 


1 


'•^f iH+,Fl.-7,6 

Fur MAT ( //^x,* ...==(= 

:3iGMA>, Fi j.6,2{iH+,Fio.6) 
’■= 0.0 

fij i-.6. KICK 
S1GMA = T/TJJ 
C0SA = CnS{SIGi'lA*Y{5)+WA ) 
COSd=:C-jS( CIGMA«Y(5 )4-W6 ) 
C‘JSC=CnS( ^-IGMASisY ( 5) + wC ) 

COS P =cn S ( J I GMA ^ Y ( 5 ) + WA +WA A ) 
CUSR=C0S{ iIGMA*Y(5 )+WB+WAA) 

tKPR^i'iPi SIGMA^YCI) ) 


'^^^2y^‘3},Y(7)!*WBD!^wiAD^BBY(5) 
.■•6f BXP ( S I GHA* j Fi 0 . 6 f * ) CCS ( 
jFi0.6tX<SIGMA+*,Fi0.6) 

SIGMA* I Fi.. . 6 j * ) CCS ( 

3.6,*£XP(SIGMA*,Fi0.6,*)CCS{ 

f r i U * D / 



t OS AD = C CS ( S IG^‘|A*Y( 7 ) +ut An » 

cjsp'Sci? 5 

C- iS R i-'-r rc ^ ^ ^ > +W AO + WAAD ) 

«ppn: 

PSI=eKP»i;.'p?tcuSptB?l5)’»s!rJJ5***‘*’**'‘*‘>*EXPRD«COSAO 
5iTI«CKCf XPR«CnSr+rci 51*|!«r^»“? 
Pul.V-:LMKi..*PvpR^o;^'^’;'r«D* = XPRo.COSCD 
Y«RVEL>PK-h<.7XPR.C isLrpI 


/PI 

. 8 0 • 

' WV-L*^S 0 . /PI 


.■Jt 


5 0 

OY; 


950 


°HID(k.I CK )=PHI ^^-0 '/PI 
PSlD(KICK)=PSl4=IrQ;/pj 

K IC.y:, ) =Hr TA=S^ -P-. 

H.JVrLDr)<lLK)=:RGLVcL«ier../pT 
YAY-LD( KIC.K)=Y ‘« — /PI 

I *( KICK )= i 
T-T+ . =c 

co'-rriGui 
'7>iJ To 6 £ 

PKINT jgo.j 

FJR*J& I ( //. 

.*46 {►) = =!. j - 

■' H’ ■- 


f UR Ma T { /5 .; (“ 5 X I" F P ^ H ^ ^ ^ ^ » A ^ < 5 } , A A { 6 ) 


c-RInT a 
FOR lATI //, 




• P t i X ) , F3 , 5 , =«=+- 1*. FR . 5 . > V ./cv r •■ c ^ 

» ro. 3 , ^ ,x, ^ { 5 X, Fd. 5 , 5 X ) ) 

= »i7X,*BB(p) jBB{4)=«,17x» 


P«jr 2 Y.O.BBIi), 3 B(EI,RBD, 81 BD.BB( 5 ),Be( 6 ) 

.*CC{5)=:*) " TA- <» ^^f-'')=*»i7X,*CC{d) ,CC(4): 

PRi'yI^ 

FdRYiAT ( // ,uX i*aA Y ( 1 ) =❖ , ; 1 X .*A A y f -. \ * 

SJ';; 

t OR MAT < // - pX, * 83 Y{ 1 ) =={Y ,i] x,)^BBY( ?)=: 

3 S ,j,b 8 Y( 1 ), 3 BY{ 2 ) jBBYIS)" 

30^7 


-ITX , 


PRINT 
PRINT 
HRINT.- 
kRInT 
PRINT 
PRI n 
PRINT 
PRINT 
PRIi^T 


'» liX,*BBY( 5 )=^) 


- - .'3 

3001 

liisiliiiiiiiiiiieSiiiiiiiiSsiji::-: 


PRINT AOu, { AAY( 1 Y / I T-' - r . ’ ' ' JtWC0,CC(5) 

print 4 - f { BS Y { I ) I Y ( I ) ,’ I = 1 7 ) J akApD ’ Y { 3 ! ' ^ ^ J ’ » WAAD , A A Y{ 5 ) 


.-) + #, Fit. 3 ,=i=SIGMA+*, FiO .6 ) » r-u.o, 

FURMAT { //;;X, ❖. 

.^XfFiO.dt 


iH+,Fi 0.6 







cr V 
- $ - 


FORM&TI _ 

■»X| rj 0.6, *tXP{ sigmU’c 
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* i f 


KILL- 

SIGM^=:T/TuU 
SUH=U, c 

SUMp = : , r 

SriP‘ , r 

* UT = 0 « 0 

fX'j i._,- 3 u { 

“ = S I GHA { 1 ) 

n=Exp(P) 

' i ) = 44( i )*Zu( I ) 


»2 


P A ,h 


YJl ) =A£Y{ I )=!cZD( I ) 


- 1>( n=B e{ I )«zo{ I ) 
ABOY(n=ZBY(n*zn{i) 
■^-{i)=cc(i)4cza(i) 

= ( I ) 

iU'^P = i(J /irD + [,;^ 4 Y ( J ) 

^ <10 = 5 ^0 +Arj { 1 ) 

SMJP=s^'IP+mGBY{ I ) 

r jT=T 0 T 4 H^^(I ) 

J ij T f :''jU E 

f u S A 0 = c c S { S I G A Y ( 7 ) + w A D ) 

.-JicD-C C i { GI GM4«Y { 7 ) +urn ) 

^ ^ I LL ) = PSI HO, / PI 
-clAC CK ILi_)=8GrA^i3 !,/pT 

J?,';-,‘-?!7“-‘-'=«aLV5L»160./PI 

11 uuu- '‘‘‘ ' 

r = T+:,. 5 

CuniL^ui 

CDNTPiiUE 
print ISO-^ 

PRINT 
PRI NT 


3 C ’ 


iJ— ) 

PRI;,T ic^; 

-05 2 i = i,_ 2 ^, 

"'-05 i J=l,:; 

TP{ J)=120 
TIMfcf j, I)=T.( 1 ) 
PARAMU. ,I 5=PHID( I ) 


<n(J,,PH,ou,,PSIO<J,,B.TAD,J,.«ovFU„l,,VAVFLO,J,, 





4^, ^ J 

2 3 C? X 
25ai 

-’00 

E365 

300i 


->‘'-i i^iTDU 


*i o n i -* n 

» A-i 0, liu, 0 , i2 , 9 ) 


»13C,120tO,i2,9) 


^ ^C, 12 : » ■■;, i. 2, 9) 


PAR a? 4( 2 ,T )=pr j j 
) = Er. T4D(l ) 

^ARAM(i.,i )=H(jv,-i_D( T ) 

PaR4M{ 5 ,1 )=YAVCLn{ I ) 

co^gTiNui 

=RI0T icX' 

PRINT 

print 3r_- 

PRINT 3C27 
CALL GR AF { Ti , d;: 7 .^3 
PRINT 3'- 7 

print 3C . 

PRINT 3-.-- '2 
PRINT 3Ci^,.3TDO 
PRINT 2 * 

PRINT 3C2/‘ 

NALL ORAF{Ti,PHM 
^RInT 3C2. 

PRINT K 
print y 
print .-.Clj,OTuQ 
PRINT 3Cu 
print 2.,2V 
Call grafj T i , psid 
print -C’, 

PRINT ICO, 

print 

print FCi ,5TD3 
print -CO, 

PRINT 

PRINT CO 
PRINT 3C_; 

PRINT 3C- ,bTDO 
PRINT 3C 
PRINT .;C2 

CALL GRAFITijYAYLLD, j.20,i2 / , . f ‘2,9) 
print iC2T 
PRINT 3Lj . 
print rC.A 
PRINT ill jbTOn 
PRINT It 
print 1C2; 

CONTINUE 
format (5313. 6) 

format ^5F7.-^) 

format ( 3F1.5) 
format { lX,i30{lH«)//) 

FORMAT (//;oA,TA=^,Fii,5,5X,*B=#,Fi 
,^D=^,F1 l.“ ,3X,#=F==i',Fii. .5//) 


,5,5X,*C=*,F11.5,5X, 




^6 


» i ' and GtOMETRICAL CHAR 

INeRTIA»i,.X*SPAN 

-•^-... A^WI NG AR-A*/-., X*(LBS) >!'i4X*( SLUGS 


H30‘ 


j46X*iXX*l 
FORMAT ( // 
:TY^iOX* 
;GS/CU.F T) ■ 
^2u.6,F17., 
FORMAT { // 


X^^'IZZ^ 


;ox 

LIFT 
■^oX^CO- 


^FL 1 uH .^^CONHIT I QNS^/P.OXi.? ( IH- ) / /30X={'V£ LOCI TY*10X*D6NS 1 
** 0X* ANGLL UF *iUX*FLlGHT PaTH^p/SZX^ { MPH ) # 6X^(SLU 
;PF^t^^-f^T*7X*ATTACK(DcG)4i .x=«'ANGLE( D£;G)«/3 2X,F6.i,F 
^jF i.o.ZtFlS.O ) 

.,.uLL R-.TL-.iX^YAW RA7t.=!=i3X^CONTROLY/3jX^CYB=«F7.4, iiX»!^C YP = *FT .4 , . 0 
.-+, i„X*CYC = ^F7.H-/3 Jx«CLB = *F7.4, liX*CLP=*F7.4,i jX=«CLR=*F7. 

"CrjB = >!'F7.4,ilX«CNP=*F7.4,l-X*CNR=4F7.4,luX*CNC = 


SC.FT)*11X*(FT)«11X*(SG.FT)*/ 


:X=^CYR = *F7 
='4,1..'<^'CLC^ 


=*‘F7.4/3 X = 



^*F7 .-) 

3007 

f^ORMAT { // 
i.ATFS(DeG/ 
2.i/29X*SI 

3x003 

FUR HAT ( // 
^40X52 ( IF- 
21*1 3 X * P £ I 
'N''*(OEG) *L. 

:005 

FORMAT ( IH 

300q 

format { IH 

0'^v9 

FORMAT ( 39 

3 " 1, 

FORMAT( 39 

3011 

FjRr'!AT{ 39 

3 ' A, 4£w 

FOR YAr { 29 

3 '33 

FORMAT { 39: 


FuRHAT( 39 

i022 

FORMAT ( 39 


FORMAT ( 29 

3027 

F0R'1AT{ 3 9. 

STOP 

;:ND 


'-rti'/T^^ COOniTlDNS*/ZuXi.d ( iH- ) //i .X^ANGL-S { DfcG ) ^iOX^cR 
‘X^RGL L*F8 . 1 , i JX*R0LL*F6 . l/idX^YAW^Fg . I t 1€X*YAW«F7 


:>iSLIP«F5. j, ) 

OX=^'VARI ATICif] 
)/54X*AHD Yaw 
«= 1 •XYHSTA^'^ •X’i'ROLL Yt 
.^(DcG) 44X*(DEG) *i2X> 


OF RULLf YAW AND SIDESLIP ANGLcS ARU ROLL*/ 
RATES WITH T IML*/54XZ3 ( IH- ) /i8X*T I RF*iwX*PH 
LGCITY*! ..X*YAW VELOCITY*/ 4.SX*(ScC) *8 
(Dtb/StC )*14X*( DEG/S FC)*) 


iH*, 


i7X,*RLjLL ANGL^ VS T1 HE* , 17X , iH* ) 
lH*,9X,*aiST'JRDANCt — BETA=*, F6 , 2 , *D£GR ELS*, ilX,iH*) 
H*,17X,*YAW AWGL'C VS TIME *,17X,iH*) 

H*, i5X,*SI DESLIP angle VS T I Ml*, iSX , 1 H* ) 
H*,14Xt*R0LLlNG VELOCITY VS T I ME* , 1 4X , iH* ) 


■ , iH*,l4X,*YAWING 
;,iH*,6X*C0EFFICl 


• , i H* 
, 54 ( 


f i6X, 

iH*) ) 


VLLOCIuY vs TIME *,i4X,iH*) 
iNTS OF CHARACTERISTIC EQUATION *6XiH*) 


*TwUATIONS OF MOTION *,i6X,lH*) 


SUB ROUT IN:: QUART (A,Y) 

THIS PRCG-’AMME called SUBROUTINE QUART SOLVES QUARTlC EQUATIONS 
BY CALL IN . TH- CUBIC AND QUADRATIC SUBROUTINES 
01 Me NS I G ; A { 5 ) , Y ( E ) , AC ( 4 ) , Au ( 3 ) , BQ ( 3 ) , RT ( 3 ) 
lQUIVALENCE( AQ tBQ) 

C3=A(3) /A{2) 

C2=a{4) /a (2) 

C]=A{5 ) /A{2) 

C„. = A(6) /A1 ;) 

u — C , j / 4: • 

aC( i ) = 1 i 

AC{/:)=-c; 

iCI 3 )=C 1*C3-4.*C ' 

AC { n ) =C C* { A . *C2-C;: **2 ) -Cl **2 

Call C'jeic[ac,rt,rti ) 

IF(RTI) e.,ii,2- 

IFtRTQ )-r;T( 2) ) il,12,:2 

Rr{I)=R7(,;) 

iF{RT{; )-RT(3) ) 13,2C>,2 . 

P.T{ i )=RT{ ) 

S=RT{i ) tZ, 



S'? 
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C 

/*> 


iOOij 

-r -j'j 5 

iOO:j 

x0i4 

^016 
:£» ^5 

1020 
’ -. -"j -% 

i ^ 

1-024' 

l'^26 

\ |> 

4; i« Cf 

i, 0 JiS ij 

1032 


IF{B=^^2-C ) 22t2£,24 

r% -rTf 
l^ ** 4f m 

CA = SQRT (C*C + 2.*R-C^) 

GO TO 25 
D=SyRT( B« 3-C , ) 

CA=- (Cl /£ .-C=^B ) /D 
2Q( I ) = : , 

&Q(2)=C-C.-. 

AQ { 3 ) =B-D 

oALL 0040 (AU,Y( 1 )jY( 2 )jY( 5 )) 

eQ( £)=C4C.- 
8Q( 2 )=34D 

BALL QU £D ( BQ » Y { 3 ) j Y ( 4 ) , Y ( 7 ) ) 
Y{5)=-Y (5 ) 

YCS )=-Y (7 ) 

PRI ^JT 6? j Y ( A ) , Y ( 3 ) , Y ( 2 ) , Y { 6 ) 

PRIM 63,Y(3),Y{7),Y{A),Y{.4) 
FORMAT {//. OX, ^THt ROOTS ARE*, Ei 5.8 
2 _ ^ 5 . S ) 

RETURN 

“ND 


2H4l,Ei5,^//2iKX,fci5.8,2H + I, 


SUB. ROUT IN . 

THIS SLB'.auTiUr- 
ROUTL.r: UjARTIC 
OIMEJSICN A(-), 
iPATH=£ 


CUBIC {A,X.R,XI) 

suLVcS Cubic tyuATioNS and is CALLtO by the sub 
and I,-i its turn calls SUBROUTINE QUAD. 
XR{7),Aa{.) 


2X=i.' /3. 

IFj .t.(4) ) I. ;6,I . :4,1 06 
XR( i )=0 .0 
CO TO 
A2= A ( 1 } **: 


y-{ i7.*42*A {4)-9.*A ( .. ) *A( 2 ) *A{ 3) 42.*A (2 )** - ) / { 54.*AZ*A( ) ) 
IF{Q) IvlO, j,00e,10..4 
Z=0.0 
GO TO iC 3 


Q 


y *“ W 

IPATH=1 

P={5.*A(i )*A(3)-A{2)**2)/(9.*A2) 

ARG = P** 24-u*Q 

IFIaRG) i-*6,i2.3,.,2- 

Z=-2.*5CRT{-P )*CnS { ATAN( SQRT(-ARG ) /Q) /3. ) 
GO TO LC2; 

Z=-2.*g**.:X 
GO TO 1C2,; 

SARG=SQRT (ARC) 

IF(P) 1. _2,1_24,av26 

Z=- {Q4SARo)**cX-(Q-S^RG)**EX 

GU TO _C2-- 

Z=-( 2.40 **rX 

GO TO 

Z={ SARG-W )44lX-(SARG+ 0 34#t:X 
GQ TO { I0_ 0,j.U22) , ICATH 
2=-Z 

XRI i ) = { 2.*A(„ )4Z-A(2) ) /(3,*A(^ ) ) 



!\.i r-o 


i.d34 aQ(1 )=&{!) 

AQ { 2 ) =A {i } +XR ( X ) *4 ( } 

,4Q { 3 ) =A (3 ) +XR ( i ) *4Q { 2 ) 

LhL L C'b 30 ( AQ I xR ( 2 ) , XR ( 3 ) . X 1 ) 
RlTURKi ^ ^ 

£ND 


:-U 

c 


30 

4 ^; 


SU3R0UTIUZ gUAD(A,XRi ,XR2,XI ) 

QUADRATIC EQUATIONS A,MD 
tH- «UAR,iC And cubic SURROUTIN-S 
D1 MwMS ! Ci j A ( 3 ) 

' ' ~ ~ “ I <- ) / I <i. • ^ A ( X ) ) 

DISC = X^ (3 ) /A( 'I ) 

IF (DISC) 

X2=S0KT {-r-iso 
aR1=X. 

XR2=x::, 

X1=X2 
GO TO 3C 

— iORT (DISC ) 

XRi =Xi+X2 
XR?=Xi-X^ 

Xi= . 

RiTURu 

,00 

SUB ROUT 10 , arctAN (X,Y,£) 

PROGRAMM- TO 3VALUATC ARCTAU Z AOD TO DET£HMIn£ THt 
PI - Y. ^ A TAO { i . ) 

Z= AT■A^ (ABS{Y/X)) 

Zu=Z*i iC. /?i 
IF ( Y ) A J } 0 } 30 
■' i F ( A ) _ 5 , 2 5 » 0 5 

5 Z=Z+P1 

T-Vj — L L' + IS.', 

GO TO c 

5 Z=2.«pi-7 

zd=36 ; . -zl; 

GO TO GC 

IF ( X ) u f ‘r 0 f 6 0 
2=PI-Z 


5 0 


ZD=„ 

s 

:■ .-z 

PR In 

T 


FORM 

-4 

T { SX 

rltur 

»-o 

5^ 

the 


PLCT 


: THt PLCT SUSRQUTPiES 'GRAF' AOD GRAPhM ARC GIV? 


IS CALLED BY BOTH 


PRUPR quadrant 


10 appendix »fc'. 
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PLOW C.HPif<T FCii; SCL,u < | O iy CF-cAT£i?Al, L-^sOAT!0<VS Of 
rA 0 i ( c N 6>v ro'.- K T H 0(< V.’ i: 1-5. Ro GC - KcTT A M £ IH BD 




fVcK- -- CM vmF (\i bW'.RALJ iC u f {MvftM tT('i?S 
'Xo 




k ^ 


K- 


t" 

k 


•6' b 

I 


+ S'-n') 

Kxz=(K2_- kL) Sx«n CcoVj 


K.- L* 




K" = Kxz 

U 2l' 
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|. 

r% vr,/^_.cy 

^/-C I 


/ 

rKiW- 

"'1-^ 

Mr. t ■ 

A-- ,c 

/" “ — 

4-^. '4. 

f) . /- 

L/ I 


* . I i. ^ j 


b. .r,. Q 

iK‘ ‘t 

,, /*-<- 
>b--, i — , r 


dc^'i 
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^ . 

bCfiSNL. 1 

" W F 

-yr T t > t '*’ ? p £ t'V ' 

"Di- 

\Xi ::: 

l-hl'2. 

i 

\ 

X^f,- 

*^'d'|>)/ Dt^f 

c. 

i- 

dn- 

'■^■^i'‘''t!.)/ /.jT fvj 


= Cl/= 

" ’'d/ I>l!'\‘ 

A 

r r t . 

■■ c) / O'i N 

r>^ 


/ I f: iV 



■ l^-4.d^')/DE.W 

■ '-3 

- 

~ Ft L^)/ ‘>):!\/ 


.(n,. 


! ^1' 

^ h 

! 

i 

- 

1 1 

' \ 

F. 


\ 

{ — o 

- C/ /.. 

\ 

E+ 

_ V, 

- <JC. 

1 

i 

f- 

-- Q’5 

■ ' i 

i t, 

i . 

{ 

= 0*S 
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N f '■< \ i At. I.NTf r ft A T I V' fv 
) Nil ) AL VALL li-S ^ 

I N \ e t?V At 

O-Off/r 
XLAlT- t-oy-f- 
X =r 

V - ^-o 
Y-'o 

*■■■ -o 



li;lJ 


. fe. '?i-) 


-V- - -*■ -J 

hZ-luHl PAT i-' AAl<,A-j 

KvCTCs^':* I 

T/ M £ii)= xr 

0(X)^ y 

^-(t) r 2 

= W 

^ (Xj = U 
AiX)^ V/r 


JZ.. 


A iJ,-r H ^ p- yii) 

A V, - Pj V Fj, / \j yy^ z) 

W , - )-'i X Fj (. U, V, X Y,' 2 j! 


lA, v: 


ki\ 


z_. 


H X V 
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j 

1 

: AU2~ ! FiO-J-t- , v-t- A\, . w v ^ vNi^.x-^ ''/i, y-+^V,, 

I AVi- ^ hJ,^ ^ x^\A!, ,'y^ >:sy^,z-i/:\zj 

/A W^-. l-lv ^ X + M/. y-, ,\y z\J:M,) 

^ Y; ' (aj-*- Avj, ) 

' ZA JZ. 7 . ^ >^ i' y \ ^ ^ J 

A Ly - 1-1 V A (g + A ij ^ ^ X/+. /^v^ _ + A Wz. ^x-+ H /z y 4 ,3 ^ z + az^ 

A V^ ~ H >< A ( U 4 - ;avz . iLJ.\y Wj , x+ hi- //+• avg Azy) 

A '^ 3 ,r FjCU-v- dUz. V-i- iivV;: ,W 4 . dWj x-t 

Y 3 “ lAx (AJ a 
^3 ^ 14 X ( V -+ ./X V 2 X 

!-};< Fd V+ZiV^ ,,V 4 - AWj X+-\Yy-f/iY 2:.-+.d.>-^ 

A Y4: U .X K,(l /v.y V- A.Vj , ■A/^/\w^_ X + w, \-(.d\y^,-~'-+ 

A 'W4 - H VA ( o -V L% gy V - V, W ^ AW^ , A -t- Vi A-" A Y3 . 

VzL- H X (u+ AU^) 

! f_ ^ ^ _ . 


X = X -t H 

y ^ Y+ i (Ay,+ 2 Ay 2 _+ 2 A Ay ) 

41- X-+ ■c:(Ai:,4 2Az:^+2A-2;3-+/:/z^ ) 

U = U-4 - i ^ 2 A v.y-f 2 AV:^ -t A ^ 4 ) 

V- V^ ^(AV,4.1A V.tIAV 3 - A V 4 .} 

w X U\vx/i -^2 A vAf a A W 3 -+.-AHl ) 
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.. r- . 


C-GRA^MS SOLVES T 
¥ FOURTH OaCiR « 
' Of ROLL, A‘4 

■Ai«^/Ci,C?,C3 ,C-^ 
A«2/o^sOi,r: ,D4 
A‘''3/€r,0 3t'3,-4, 
^ TiMEf 5,135) , PA 
^4PH10( iJ :■) ,?3iCl, 


GPHM/iSCAt 


: ■ v L t Wk: I T ,X I KO , ZI K3, SPAN , AREA 
!f 31 1 PH I G , DPH 1 0 , P S ! C , DP S i n , P. in 
r {5F7.^) 

; 1*1,2 

; T C I , VE L M , RHO , S P S C , A L F A D , 3AMMA 
(SFlG.o) 

-LdK=l,2 

^5€i,CL8,CHB,CY3 
5t L,CLP,CMP,Cyp 
!5.C1,CLR,CNR,CYR 
. . f,CLC,CNc,Gyc 
■4 3Fd.5) 

t :x,i3oaH^)// ) 


./PI. 

;&«^£IT/5RTY 
(Al r4D~EPS0)/PID 
i4*3#Jrf»iAO/PID 
i-=?TC.J7PID 
^?St>*3S./60. 

. - ^ : T/ { A ) 

.• "''.T'.^TAMCGAWMA) 

■■.• ' :?X1K0/AMA5S ) 

. " ZI^U/AMASS) 

■' 1 '< /^CCStcTA) )^=^^4{ZK0-5I'H i/4)}=^=4.)/{ 

^ .^COSIETA) + l <Ki-4=Si.a T^) )-*.)/( 
!• ;■•' **2-Xi<0**2)*SlH{tT-^3 4'CjS( £Ta) )/{S-= 
: / <K 
' ’■•■ / IK 

f /fRHa=^AREA^Y£L) 

;• /CRHG^AR5i*S>>& ;) 

. > HA 't/XXK^CLS 
CL^/XXK 
CLR/XXK 
' - G/XXK»CLC 
■ - y/ZZK^CN8 
. : Cm/llK 

tm/iiK 
U/ZZK^CflC 


f'<j *■£ 

A.X'i PP.£-[i 
ANj RuLl 


y - L L { 2 
I z 


S I' .4 N * ♦ i. ) 
S P i iN’ * ^ 2 j 
A » ) 


Cf7 


^ , 'Fa,. -• i 3- 

1 •' T > 1 h C A > >: - 

.-■A Ylv ■\-M = :r. 


) , Y A ^ „ L C ( ^ ) , 





^«iTAQ 











V+DELV3 , W+0-LW3 
V+0ELV3 f W+DELWi 
^ + OELVE, W+Oc LW3 


f X + H, Y + Ot LY? »Z + DELZ3 
fX + H,Y+D£L*Y3,Z + DELZi 

f X+U , Y+Ul LY- f Z + CZLZi 
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table 1 ; MASS GHARACJEfilSTiaS OF OESSH ArlggH 


■M 


aramet er 


Units 
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\~RT.E 2 : VARIATION 0? FUSELAGE CROSS-SlpTIOHAI) 
AREA. ALORG BODY OMTRE LIN6 


on Distanc© from Cross—Secf ion Avorago Gross* 

Station Zero (in. ) Area (Sq.Et. ) (b)xCc) sectional 

Ar ea~ 




9 . 5 

4.58 

43.5 

39.0 

13.78 

537.2 

69.0 

18.20 

1266.0 

105.0 

16.92 

1779.0 

132.0 

13.00 

1716.0 

150.0 

9.44 

1416.6 

213.0 

4.33 

922.2 

253.0 

1.10 

278.5 

=970.5 

z(d)=7960.5 


y''4x(e) 

verage diameter d^^ - ^ 


=3.26 ft, 


8.3 Sq;rt. 
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- POF? CESSNA- 182H 


Component Parameter 

Value 

Units 

liM Span b 

3 6. 083 

ft. 

Area S 

174.38 

Sq.ft. 

Aspect Ratio A 

7.46 

••• 

Taper Ratio ?[ 

0,672 

.... 

Section C^ 

, a 

0.105 

(deg""^) 

Sweep back angle 


0 

deg. 

M.A.C, c 

4.925 

ft. 

Position of C,G. 

26.4% of 



M.A.C. 


Dihedral I 

1.75 

deg. 

Horizontal Snan b.. 

11 .667 

ft. 

Tail 

Area Sp 

41.3 

Sq.ft, 

Aspect Ratio Ajj 

3.33 


- . Taper Ratio/^'jj 

0.59 


Gi (Section) 

0.105 

(deg"”^) 

a 

A: A , A' 

! c/2’ c/4 

0 

deg 


3.667 

ft. 

Th - 

0 

deg 


Remarks 
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Parameter 


Span b-y 
iin 0a. 

Aspect Ratio 


%/| 




Value 


5.333 
19.36 
1 .468 

0.526 


-0.862 

24 . 1 25 
165.0 
0.57 

0.69 


Unit 


ft. 

Sq.ft. 



n 

a 

0.105 

(degT^ ) 


/\ 

45 

deg. 


-\/A 

40 

deg. 


^/2 

36 

tio 

CD 

nzS 


> 

1 o 

3 . 667 

ft. 



2.4166 

ft. 

Body 


-2.083 

ft. 


-3 

4.833 

ft. 


Remarks 


ft. 

in. Big. 5 

mm 

- Pig. 6 
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(G0g3)I)« ) 


Unit 


Remarks 


?alue 


Sq.ft. Wei ght ed 

average from 
fable 2. 
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TABLE 3 ; (QONG?D . ) 


Gomponent Parameter 

Value 

Unit Remarks 

Other Pactors 

1.72 

- 

Pig *4 

k 

1.0 

- 

Pig. 7" 

tusi. 

1.2 

— 

Pig, 8 





fvdsl 

A^(B) 

1.0 

- 

Pig.9 

1.14 


Pig, 10 


0.0 

(deg.“ 

’’) Pig, 11 


0,0 

(deg.- 

Pig, 11 

V 

1.0 

mm 

Pig, 12 

A 





0.99 

•mm 

Pig.13 

(Kf^H 

1.00 

— 

Big, 13 

((Cl 

P 

-0,0005 

(deg.“ 

'^) Pig, 14 

(( ” 

-0.0021 

('•) 

Pig. 14 

(Ci/r^w 

-0.00023 

(deg.“ 

'^) rig. 15 

P 




V ' Ai 

-0.000155 

i ( »» 

) Pig. 15 
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TABLE 3 ; (CONTD.) 


Component Parameter 

Value 

Units 

Remarks 

(% p V 

1.00 

- 

Rig. 16 


1.00 

- 

Rig. 16 

(AO-, /O tan 

P 

c/4^W -O- 000037 

{aeg-2) 

Rig. 17 

^ » n 

)jj -0.000025 

ff 

Rig. 17 

% 

.00115 

- 

Rig. 19 

31 

0.9584 

MV 


(P'A/v)^ 

7.79 

- 


(pVKJh 

3.45 

- 


P Cl 


(rad*"^ ) 


^ K 

-0.485 

Rig. 23 


-0,24 

( " ) 

Rig. 23 

(^) 

Cl 

-0.1085 

( ” ) 


11 

o 




o 

o 

-0.1085 

( " ) 



M 


1 

(ZiG /e) 0,000057 {.rad-d^X 

Up 


^ig.25 
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Table 4 : flight 

CONDITIONS iiND ALLIED P aRAMF'PT^^R.q 


Parameters 

Flight 

Oonditions 

Velocity 

MPH 

iHtitude 

Ft. 

Slugs/cu.ft. 

Approach 

73 

S.L. 

.0002378 

Cruising 

150 

5,000 

.0002049 


TABLE 5 : PAR/J\fETERS LEPENDEHT ON FLIGHT OnW nTT tawo 





XUXNQ* 

Parameters 

Unit 

Approach 

Cruising 

^L 

- 

1.1150 

0,306 

Od 

- 

0.132 

0,032 

s 

a 

- 

0.547 

0.121 • 


— 

0.3975 

0.01393 

It.'V 

ft. ^/sec • 

156.4x10“^ 

176.6X lO""^ 

B 

^e ^ ~ 

— 

16.51x10^ 

30,15x10^ 

\ 

- 

1.61 

1.74 
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TABLE 6(a) : GOaCP ARISON OP CONTRIBUTION OP VARIOUS COMPONENTS 
TO STABILITY DERIVATIVES ESTDUTED PROM 
REFERENCES 1 TO 5 (ilPPROACH aT 73 mph) 


Deriva- 

tive 

Ref.1 

Ref .2 

Ref .3 

Ref. 4 

Ref .5 

Manufacturers Data 

Pw 

N/A 

N/A 

N/A 

N/A 

-0.0099 

N/a 

4- 

Cy 

n/a 

N/A 

N/A 

N/A • 

-0.1186 

N/a 

^B 

Pv 

N/a 

N/A 

IJ/A 

N/A 

-O.I 9 O 8 

N/A 

Gy 

^6 

N/A 

N/A 

N/A 

N/a 

-0.3191 

- 0.303 


-0.1208 

-0.1208 

-0.1817 

- 0 . 01208 - 0. 0919 

N/a 

%B 


Neg. 

-0.0028 

Neg. 

Neg, 

-0.0024 

N/A 

Ph 


0.0026 

0.0024 

0,0026 

0.0024 

0,0024 

N/a 

Pv 

°1 

-0.1182 

-0.1212 

-0.1791 

-0.1183 

- 0.0919 

- 0,097 

CCL 

0.0249 

0,0274 

Neg, 

0.0249 

Neg, 

n/a 

C 

""Pb 

-0.0492 

-0.0492 

-0.0012 

- 0.0492 

-0.0276 

n/a 

G. 

0.0928 

0.0928 

0.0928 

0.0928 

0.0857 

n/a 

% 

C 

0.0684 

0.0709 

0,0916 

0,0684 

0.0580 

0,0701 
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R6f*1 Ref » 2 Ref •3 Ref »4 Ref#5 Manufacturers 

Data 


0.0004 


0,0006 Neg 


■ 0.0022 - 0.0022 


0.0028 -0,0022 Neg 




.nufacturer 


No method available or value not available 


2» The values of roll“rate derivatives estimated by the 
method of Ref, 6 are as follows; 
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Ti3LE Sfb) : GOlgi.RISON OF CONTRIBUTION OF VARIOUS 
GOI£PO HENTS TO STABILITY DERIVi.TiVES 
ESTIM ATED gROM R EFERMCJ5S 1 TO 5 
(' CRUISE AT 150 m,p.^h ) 


Iterivat ive 

Ref.l 

Rei .2 

Ref .5 

Ref .4 

Ref. 5 

Manof act- 
urer’s Data 

%W 

W/A 

R/A 

R/A 

R/A 

- 0,0099 

R/A 


N/ii 

R/^x 

R/A 

R/A 

-0.1184 

R/A 

■ ■ 

AiftiM- ■ 

N/a 

R/A 

R/A 

R/A 

-0.1908 

R/A 

i»s 

R/A 

R/;* 

R/A 

R/A 

-0.5191 

- 0,303 


-0.0884 

-0.0884 

-0,0675 

-0.0884 

-0,0688 

R/A 

■ 

:.:-5Ph 

Reg. 

0,0015 

Reg. 

. Reg, 

0.0015 

R/A ^ 

■*fl 

Py 

-0.0155 

-0.0125 

-0.0155 

-0.0125 

-0.0123 

,R/A 

: P 

-0.1017 

-0.1006 

-0.0808 

-0.1007 

-0,0798 

- 0,0923 

... G . 

■ Pw 

0.0019 

0,0026 

Reg. 

0.0019 

Reg, 

R/a ~ . 

;. Pb 

-0.0492 

-0.0492 

-0.0012 

- 0.0492 

- 0,0299 

R/A 


0,0918 

0.0918 

0.0918 

0,0918 

0,0849 

R/A 


0,0445 

0.0452 

0.0907 

0,0445 

0.0550 

0,0587 





Ref«5 Manufacturer* 
Dut a 


0,461 

-0,4596 

N/ 

0,0033 

-0,0033 

W- 

0,0001 

-0,0017 

N/, 

0,4643 

-0*4646 

■*0 II 
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No method ayailable or value not a'’ 


e 2 ::Ihe values of roll-rate deriva 
^'imethod of Ref. 6 are as follows 




Ref.l 

Ref. 2 

Ref .5 

Ref ,4 

Ref, 5 

Manufactur- 
er's Data 

=0.0755 

-0.0755 

-0.0755 

-0.0755 

-0.0755 

N/A 

-0.0864 

-0.0864 

-0.0864 

-0.0864 

-0,0864 

-0.0937 
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PMBiiaaM^ 




not estiEiat ed 

less than the manufacturer's value 
rr,nr>P thaii the manufacturer's value 


GOICP..RISON Of ESflM/^TEI) AND AVAILilBlE 
VALUBS OF SfhBIlIfY DERIVATiyES (hPPIlO.iOH 
AT 73 m.p.h.) 


Percentage Error From the 
Manufacturer's Values 



JAILABLE VAJUES 


GOtlP.JlISOH OF ESTIMATED I^D 
OE STABILITI TlERI VATIVES (CRUISE AT 


From the Manufacturer's Values 


N,E, means not estimated. 

» mean, less than manufacturer 'a value 
*♦ meana more than manufacturer's value 




BAILABLE 


A CONoOLIBii^TED LIST OE ESTIMED. 
YiUiUES OE stability BERIVATIVES 


Approach 


condition 


Est imat ed Manuf ac tu- 
ner’s Data 


Estimated 
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PERIOD DaMPlEG GHARAGPERISPIQS 


PiPlE 9(a) 


Root 


Manufac- Estima- Manufac 
turer's ted turer's 
Data Data 


Estima' 

ted 


Dutch 


WD Di'MPIIfG CHARAGTERISTIOS 
at 1^Q_mj£^h) 


PERIOD . 
(GRUISE 


Data 
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